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Abstract 

The question of finding solutions to given implicit differential equations (IDE) has 
been answered by several authors in the last few years, using different approaches, 
in an algebraic and also a geometric setting. Many of those results assume in one 
way or another that the subimmersion theorem can be applied at several stages of 
the reduction algorithm, which, roughly speaking, allows to reduce a given IDE to 
a collection of ODE depending on parameters. The main purpose of the present 
paper is to improve some of the known results by introducing at each stage of the 
reduction algorithm a desingularization of the manifolds with singularities that 
may appear when the subimmersion theorem cannot be applied. This can be done 
for analytic IDE by using some fundamental results on subanalytic subsets and 
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desingularization of closed subanalytic subsets due mainly to Lojasiewicz, Hiron- 
aka, Gabrielov, Hardt, Bierstone, Milman and Sussmann, among others. We will 
show how this approach helps to understand the dynamics given by the Lagrange- 
D'Alembert-Poincare equations for the symmetric elastic sphere. 

Keywords: Implicit differential equations; Differential- algebraic systems; Desin- 
gularization; Nonholonomic systems. 

1 Introduction 

Implicit differential equations (j){x, x) = (IDE) are very common in science and technol- 
ogy. The case of an ODE x = f{x) is the simplest particular case. The Euler-Lagrange 
equations for a given Lagrangian, degenerate or not Pfl El Ej; Lagrange- D'Alembert 
equations for a nonholonomic system and their reduced versions Lagrange-D'Alembert- 
Poincare equations [H El EI; some equations in the Dirac theory of constraints jTj; 
electrical circuits called descriptor or semistate systems jH]; electric power systems P; 
nonstandard singularly perturbed systems ^Hj; reactive columns ^T|; constrained robot 
systems are some of the examples. 

In the literature, IDE are often called differential algebraic systems (DAS) in the case 
in which the subimmersion theorem can be applied, in one way or another, at each 
stage of the so called reduction algorithm, ^H]- In the last few years this approach 
to finding the solutions to a given DAS in the category of C°° manifolds and maps, 
and even some more general IDE, has been intensively studied. Invariants like the so 
called index of the system were introduced and calculated, as far as it is possible, for 
a general DAS and geometrically motivated algorithms have been developed, see, for 
instance, |T31 CH US UEl UH UH] just to mention part of the relevant bibliography more 
or less connected to the present work. This kind of approach has now reached a nearly 
optimal answer to certain basic questions, for instance, the question about existence and 
uniqueness of solutions for DAS. 

One may say that the reduction algorithm, as exposed in the references cited above, uses 
some of the old ideas of the Dirac theory of constraints [7], but generalizes them to 
make them useful in a variety of fields beyond Hamiltonian and quantum mechanics, like 
control theory and nonholonomic mechanics, as it appears, for instance, in the references 
cited above. 

To the best of our knowledge, the singular cases where the subimmersion theorem can- 
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not be applied have not been systematically and fully studied in the existing literature, 
although there are several related papers. For instance, an algebraic approach relying 
on complex algebraic geometry is presented in |I9|. In that article the singular cases 
where the subimmersion theorem cannot be applied are treated in the case of implicit 
differential equations given by complex polynomial relations. Also an implementation 
using computer algebra systems is provided. Besides [12], cited above, there are several 
papers studying specific questions about singular points and impasse points, see, for in- 
stance, |2D1 1^ 123 12S] • See also Jlj for examples in nonholonomic mechanics where a 
simple example of desingularization appears. 

Desingularization ideas appear in different contexts in mathematics. The fundamental 
theorem of Hironaka j21I, whose proof was simplified and computationally implemented 
in subsequent works [23 120]; is a good example of a desingularization procedure, in 
this case desingularization of certain algebraic varieties. Another fundamental theorem 
on desingularization was proven by Bierstone and Milman [2Z], which actually includes 
Hironaka's theorem. See also [2H]- 

In the present work we are going to use some of the known desingularization results, 
and also several results from the theory of subanalytic sets. Our main result shows how 
to reduce a given real analytic IDE to a real analytic IDE of constant rank, defined in 
section 121 which is considered the simplest case in this paper. We do this by conveniently 
modifying the usual reduction algorithm by including a desingularization of all the sin- 
gular manifolds that may appear at each stage of the algorithm. The only main result 
on desingularization that we use for doing this is the theorem of Hironaka on desingu- 
larization of closed subanalytic sets [22]- In fact, the theorem on desingularization of 
closed analytic subsets, theorem 5.1 of Bierstone and Milman in [^, is enough for our 
purposes. We will keep the conceptual framework as basic as possible throughout the 
paper, avoiding the unnecessary usage of schemes. 

Control systems in the category of subanalytic sets have been recently studied in [^ . 
where desingularization techniques, exposed, for instance, in |30^ and [12], have been 
used. In the present work we also work, in a sense, with control systems in the category 
of subanalytic sets, but from the point of view of IDE, which in a sense is dual to the 
point of view of control theory. In fact, a control system is, roughly, a vector field de- 
pending on a parameter, or equivalently, a family of vector fields, or, more generally, a 
family of local vector fields. On the other hand, as we will see, an IDE of constant rank 
gives also a family of vector fields, but defined implicitly. 

Working with IDE in the subanalytic category rather than in some of the C^, k = 1, 2...oo 
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categories, of course is a limitation, but there are interesting examples that can be stud- 
ied in this context. In fact, in the study of IDE representing several important examples 
from mechanics, control theory and other fields, as we have mentioned before, equations 
are often given by real analytic functions. 

The relevance of the desingularization method that we propose to help finding the solu- 
tions to a given IDE depends, in part, on a sufficient knowledge of the geometry of the 
resulting desingularizing manifold. For instance, if in a given example the final system 
obtained by the methods of the present paper is simply a vector field, existence of equi- 
librium points may be directly related to the topology of the desingularizing manifold. 

We develop an interesting example showing how the desingularization method helps to 
solve a mechanical system, namely, the rolling symmetric elastic sphere. By definition, 
this system has an extra nonholonomic condition besides the usual nonsliding condition 
for a rigid rolling sphere with only one point of contact with the plane. This extra con- 
dition says that the vertical component of the angular velocity should be 0. A physical 
situation corresponding approximately to this model is that of an elastic sphere, like, for 
instance, a rubber sphere, which is slightly deformed as it rolls on the fioor. We assume 
that the deformation of the sphere takes place only on that part of the sphere in contact 
with the fioor. Then there is an area of contact, rather than a point of contact, which 
imposes, because of the nonsliding condition, the condition that the vertical component 
of the angular velocity is 0. From the point of view of reduction theory, this condition 
introduces singularities in the reduced system of equations, since the standard dimension 
assumption for the Lagrange-D'Alembert-Poincare equations is only partially satisfied, 
see |3] . The nonholonomic condition that the vertical component of the angular velocity 
is has been considered for the example of a rigid body with a fixed point in We 
must say that the model that we are considering for a rolling symmetric elastic sphere 
is not necessarily very realistic because it does not takes into account, for instance, the 
deformation that may occur in the material far from the area of contact, or the vis- 
coelasticity, which introduces internal friction. Those are delicate problems in elasticity 
theory which have been the subject of several investigations for more than a century and 
are an active area of research in our days, but this is not the subject of our paper. We 
will concentrate only on the dynamical equations for a rolling sphere under the specified 
nonholonomic conditions. This kind of system has potential applications to robotics. 
The desingularization process can be performed in detail in the case of the symmetric 
elastic sphere and the desingularizing manifolds can be identified, which leads directly 
to solving the system by quadratures. In fact, we show that the original system is equiv- 
alent to a system on the nonsingular manifold S"^ x S^. In spite of the vast literature on 
the subject of nonholonomic systems with rolling constraints we believe that our results 
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in the specific case of tlie symmetric elastic spliere are new and interesting by tliemselves. 

We will work mainly in the category of subanalytic sets and maps. Manifolds will be 
usually real analytic manifolds and maps from one manifold to another will be real an- 
alytic maps, although some of the statement are also valid in a more general context. 

In section |21 we explain some basic facts about IDE. In section we describe our algo- 
rithm. In section m we prove our main results. In section |S1 we give a detailed description 
of the example of the symmetric elastic sphere. 

2 Implicit differential equations 

Basic notation. In this section manifolds will be smooth manifolds and maps from 
one manifold to another will be smooth maps. Let M be a given manifold of dimension 
n and F a vector space of dimension m. Let a : TM — > F be a map such that, for each 
(x, x) e TM, a(a;, x) = a{x)x = a(x)(x, x) is linear in x. Let / : M — >• F be a given map. 
We will study IDE of the type 

a{x)x = f{x) (2.1) 



By introducing the trivial vector bundle M x F we can think of a as representing a 
vector bundle map 

a : TM ^ M x F 

and of / as being a section of M x F. Then, for each x G M, a{x) is a linear map 
depending smoothly on x from the tangent space T^M into the fiber {x, F) of the triv- 
ial bundle. More generally, we may consider a general vector bundle with base M, say 
71 : F —>■ M and an IDE like ()2.1|) where now a : TM — ^ F is a vector bundle map and / 
is a section of F. This kind of generalization is important to describe a sufficiently wide 
class of IDE. However, in the present paper we shall describe only the trivial bundle case, 
for simplicity and also because it already contains the essential facts. The general case 
can be treated in an essentially similar way. In this paper the manifold M is called the 
domain and the space, or more generally, vector bundle F, is called the range of the IDE. 

Given an IDE one has immediately a linear algebraic system (LAS) for each x G M, 
depending smoothly on x, where the unknown is {x, x), for each x G M. We will call it 
the LAS associated to the given IDE. 
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IDE of constant rank. Assume that the LAS associated to (|2.H) has solution {x, x) 
for each x G M. Then (12.11) defines an affine distribution, generally singular, on M. If, 
in addition, ranka(a;) = rank[a(a;), /(x)] is locally constant, that is, it is constant on 
each connected component of M, then the IDE is called an IDE of constant rank. This 
is equivalent to saying that the corresponding affine distribution has constant rank on 
each connected component of M. For instance, if m = n and a{x) is invertible for all 
X E M then (j2.H) is equivalent to an ODE 

X = a{x)-^f{x) (2.2) 

and the rank of the affine distribution is in this case. 



The case of an IDE of constant rank is the simplest case in our context and our main 
result, in sectional shows that a given analytic IDE can be reduced to a finite collection 
of analytic IDE of constant rank, which can be considered also as a single IDE of constant 
rank in an obvious way. 



Reduction of a general IDE to an IDE of the type (j2.1j) . It is easy to see that, 
from the point of view of a general theory where the dimension does not plays an essential 
role, IDE of the type 

x) = 

where the map (p : TM — ^ F may be nonlinear in x are not more general than ()2.1|) . 

In fact, let us assume first, for simplicity, that M is an open subset of a finite dimensional 
vector space E. An IDE of the type 

x) = 

can be rewritten in the form (j2.H) with domain M x E and range E x F a.s follows 



X 



u 



which has the form 1)2.11) with 

a(x, u) = 



I 




[x,u} 



f{x,u) 



u 

(x,u) 



The case of a general manifold M and a fiber preserving map (p : TM F where F 
is a vector bundle with base M can be also reduced to the form 1)2.11) by an essentially 
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similar procedure, using a geometric construction involving pull-backs of bundles. 

Remark, (i) We must remark that working with IDE written in the form ()2.H) is an 
important ingredient of our algorithm, in part because this form is preserved and the 
space F remains the same (or the vector bundle is replaced by a pull-back vector bundle) 
at each stage of the algorithm, which simplifies matters as will become evident later, 
(ii) The system (j2.H) can be written equivalently as follows 

a{x)x = tf{x) 
i=l, 

where t = t{s). The first equation of this system can be written equivalently as follows 

b{y)y = (2.3) 
where y = {x,t), h{y) E L{T{M x M), F), h{y) = [a(x), -f{x)]. 

Of course ()2.3|) is an IDE whose associated LAS has a solution for each y, but it does not 
seem that questions like reachability for (|2.1|) could be easily reduced to easily solvable 
corresponding questions for (j2.3|) . In other words, this kind of transformation of the 
system does not necessarily really simplifies the problems related to a given IDE. On the 
other hand, systems like ()2.3p are interesting by themselves and are related to Pfaffian 
systems jSl] . 

Some notation and operations with IDE. It will be convenient to denote (a, /) 
the IDE (j2.H) . from now on. Let (a, /) be a given IDE, say 

a:TM ^ F, f : M ^ F. 

Let C M be a given submanifold. The restriction {a,f)\N, also written {a\N, f\N), 
of (a,/) to is defined naturally by the conditions {a\N){x){x,x) = a{x){x,x) and 
{f\N){x) = /(x), for all {x,x) E TN. 

Let us remark that the notion of restriction (a, /)|A^ makes sense also in the case where 
A^ is any subset of M. In fact, we only need to give a meaning to the notion of a tangent 
vector {xo,x) at a point Xq E N. It is simply the tangent vector in T^^M to a smooth 
curve x(t), t E (—5, 5) in M such that x{t) E N, for all t E (—5, 6) and x(0) = Xq. 

Let (f : N —>■ M he a. given map. Then the puU-back ip*{a,f) = {ip*a,ip*f) is the 
IDE with domain A^ and range F defined by {p*a{y){y,y) = a{{p{y)) {T{p{y,y)) and 
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li g : F ^ G is a linear map we define the projection of (a, /) by g as being tlie 
IDE witli domain M and range G defined by {g o a, g o f). More precisely, we have 
{g o a){x){x, x) = g{a{x){x, x)) and {g o f){x) = g{f{x)), for all {x, x) G TM. 

One can define operations like the direct sum © or tensor product of IDE in a nat- 
ural way. For instance, if (a^, fi), i = 1,2 are given IDE with domain M and range Fi, 
i = 1,2, then we can define the direct sum (ai,/i) © (02, /2) = (oi © 0.2, fi © /2) as 
an IDE with domain M and range Fi © F2 by (ai © a2){x)x = ai{x)x © a2{x)x, and 
(/i © (x) = fi{x) © f2{x), for all (x, x) G TM. The tensor product is also defined in 
a natural way. 

Using operations like the ones described above, one may sometimes simplify a given IDE. 
For instance, working in coordinates in F, say {yi,y2, ■■■,ym), if some of the equations, 
say corresponding to yi, is a linear consequence of the others then it can be eliminated 
by using a projection g{yi,y2, ■■■,ym) = {y2, ■■■,ym), and the resulting system will be 
equivalent to the given one. We will need the following result, whose proof is not difficult. 

Theorem 2.1 Let {a, /) be a given IDE with domain M and range F and let N C M be 
a given submanifold defined regularly by equations = 0, where (p : M H and H is a 
finite dimensional vector space. Then the restriction {a, f)\N has the same solutions as 
the IDE (a ©0, /©</)) with domain M and range F®H. It also has the same solutions as 
the IDE {a®DipQ)0, fQ)0®ip) with domain M and range F®H®H. Here Dip : TM H 
is defined by Dip = p2° Tip, where p2 : H x H ^ H is the projection on the second factor 
and TH = H X H. 

This theorem is simple but useful. For instance, it allows sometimes to replace a given 
IDE by an equivalent IDE whose domain and range are vector spaces, which sometimes 
simplifies practical calculations avoiding the usage of local charts whenever it is conve- 
nient. More precisely, let a given IDE (a, /) having domain M C L, imbedded in the 
vector space L and defined regularly by an equation = 0, where ip : L ^ H, and let a 
be defined by a restriction a = A\TM, where A : TL F. Then, according to theorem 
12. H one can work equivalently with the system (A © 0, / © whose domain and range 
are vector spaces. 

One can obviously define a category whose objects are of the type (M, F, (a, /)) , where 
M is a manifold, F is a vector bundle over M and (a, /) is an IDE with domain M and 
range F. A morphism ip : (M, F, (a, /)) —>■ {N, G, (6, g)) is given by a map pd '■ M ^ N 
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and a vector bundle map ip^. : F ^ G over Lp^ such that, for any (x, x) G M, we have 
iPr {a{x)x) = b{y)y, and g{y) = ipr if{x)) , where {y,y) = Tipd{x,x). However, although 
this perspective is interesting, we will not use the categorical context in this paper since, 
as we have said before, we want to keep the context as basic as possible. 

Basic reduction algorithm for solving an IDE. We have the following basic algo- 
rithm to solve ()2.H) recursively. It is a reformulation in our context of essentially known 
ideas contained in the references given before. However, we must remark that the fact 
that we have chosen to write a given IDE in the form (j2.ip where the space (or more 
generally, vector bundle) F, remains unchanged (or is replaced by a pull-back bundle) 
throughout the reduction process, briefly described next, presents some clear technical 
advantages. 

(01) Ml C M is the subset of all x G M such that the LAS (j2.ip has solution. 

(02) For k = 1,2, .. . we assume that Mk is a submanifold and Mk+i C Mk is the 
subset of all x G Mk such that the LAS ()2.H) has a solution (x, x) G TM^. 

If in a given example the assumptions that M^ is a submanifold made at each stage of the 
previous algorithm are satisfied then the algorithm itself stabilizes at a certain stage q, 
that is Mg = Mg+i. Then the system 1)2.11) restricted to Mg has a solution {x,x) G TMg, 
as a LAS, for each x G Mg. Thus we obtain an IDE with domain Mg and range F which 
may or may not be of constant rank according to the definition given before, but it is 
equivalent to the original one and the associated LAS has solution for all x G Mg. The 
assumption that M^+i C M^ is a submanifold is obtained in practice by using, somehow, 
the subimmersion theorem or equivalent results. An important case in which this scheme 
works is the case where M is a vector space, a{x) is independent of x and f{x) is linear 
in X. This kind of situation was studied in [33]. It can be also studied with our methods, 
leading to perfectly identifiable and simple desingularizing systems, which we are not 
going to explain in the present work. 



3 Desingularization 

The hypothesis that M^ is a submanifold at each stage of the basic algorithm described 
above is too restrictive as it is not satisfied for many examples of interest. Some rather 
general geometric approaches to IDE in the C°° category, like in [121 CH] , also assume 
some kind of geometric version of the subimmersion theorem. We are going to show 
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that in order to overcome part of those hmitations, and at the cost of working in the 
subanalytic category rather than the C°° category, one can use results from real ana- 
lytic desingularization theory. This is possible thanks to the now well established theory 
of semianalytic and subanaliytic sets developed originally by Lojasiewicz [321 EH EH]- 
Important results in this field and systematic expositions using techniques which are 
simpler than the original ones are due to Gabrielov , Hironaka [21] , Hardt |4U|, 1^ , 
Bierstone and Milman [30J, Sussmann [S2], and others. Our main reference will be 
|30j . where an excellent and very readable exposition of important points of the theory 
of subanalytic sets has been written. 

For the rest of this paper manifolds and maps will be real analytic, unless otherwise 
specified. For instance, if (a, /) is a given IDE with domain M then M will be a real 
analytic manifold and a, / will be real analytic maps. 

Definition 3.1 Let M he a real analytic manifold and let X he a closed suhanalytic 
suhset of M. A desingularization of X is a a real proper analytic map f : N ^ M such 
that f{N) = X, where N is a real analytic manifold of the same dimension as X. 

This is a relatively weak notion of desingularization, but it is enough for our purposes. 
Existence of desingularizations f : N ^ M is guaranteed by the following theorem of 
Hironaka, see |301 |H2] and references therein. 

Theorem 3.2 Let M he a real analytic manifold and letX he a closed suhanalytic suhset. 
Then there is a desingularization f : N ^ M of X. 

Desingularization results that include those of Hironaka have been recently proved in 
Bierstone and Milman [77j . 

In fact, in the present paper we only need the following weaker desingularization result, 
which is theorem 5.1 in [30^, 

Theorem 3.3 Let M he a real analytic manifold and let X he a closed analytic suhset. 
Then there is a desingularization f : N ^ M of X. 

Description of the algorithm. Let M be a manifold of dimension d and let (a, /) 

be a given IDE with domain M and range F. The basic result proved in this paper to 
solve the IDE (12. ip consists, roughly, in transforming it into an equivalent IDE, say 

a2{y)y = My) 
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on a manifold M2, which is a DAS of constant rank. The manifold M2 will be constructed 
by an algorithm that involves a desingularization process. 

The decomposition M = Mq U Mi U M2. First, let us assume that M is a connected 
manifold of dimension d. For i = 0,1, . . . , let 

Si{M) = {x e M I ranka(x) < i} 

= {x G M I det A{x) = 0, A{x) submatrix of a{x) of order i + 1} 

Si{M) is clearly a closed analytic subset of M, defined by analytic equations, for i = 
0,1,.... 

Also, for i = 0,l, let Li{M) C Si{M) be defined by 

Li{M) = {xe Si{M) I rank[a(x),/(a;)] < i} 

= {x E Si{M) I det A{x) = 0, A{x) submatrix of [a(x), f{x)] of order i + 1}. 

Each Li{M) is a closed analytic subset of M defined by analytic equations. 
Let 

SkAM)cSk,iM)c...cSk^iM), 

where kr = kr{M), be the distinct nonempty Si{M). We observe that Sk^{M) = M. 
Consider the corresponding inclusions 

i:fc,(M)CLfc,(M)C...CL,,(M). 

We have that ranka(x) = rank[a(x), /(x)] = kj for each x G L^.^M) — Sk _i{M), 
j = 1, . . . ,r. The LAS associated to ()2.1|) has solution for each x G Lk^{M) — Sk^^^i^M), 
j = 1, . . . ,r, where we have, by definition, Sk^ = 0. 

We remark the following useful facts: the set Lk {M) — Sk -i{M) may be empty, for some 
j = 1, . . . ,r; we have dim Sk,_^{M) < dimM; if dim(Lfc/(M)) = d, then Lk^{M) = M. 

Now let M be a manifold of dimension d and assume that 

j 

is the union of the connected components of M of maximal dimension d. 

We will consider the following pairwise disjoint conditions for a given Wj C Mm, 
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(a) LkXWj) = ^ 

(b) LkXW,) ^ and ^imLk^Wj) < d 

(c) LkXWj) ^ and dimLuXWj) = d. 

We now define the following pairwise disjoint subsets of M. 

Mo = (M - M^) u\jLk^{W,) u\JSk^_,{W,) 

b c 

M, = [jw,yj[j{w,- L,xw,)) 

a b 

M, = [j{W,-S,^_,{W,)). 

c 

We have the following assertions, whose proof is easy: each subset L^^iyVj) C Wj, and 
each subset Sk^_^{Wj) C Wj, is a closed analytic subset of Wj defined by analytic equa- 
tions on Wj. In consequence, Wj — Lk^(Wj), Wj — Sk^_-i^(Wj) are open submanifolds of Wj. 

The manifold M is the disjoint union 

M = Mo U Ml U M2. 

The manifolds Mi and M2 are open submanifolds of M. The subset Mq is a union of 
subsets defined by analytic equations on each Wj, union M — M^, and we have that 
dim Mo < d. 

Restrictions (a, /)|Mo, (a, /)|Mi, (a, /)|M2 and desingularization of (a, f)\Mo. 
We have that the LAS associated to ()2.H) has no solution for x G Mi. On the other hand, 
it has solution for all x G M2, moreover, (a, /)|M2, is an IDE of constant rank. 

It remains to see what happens with the system restricted to Mo. The idea here is 
to desingularize each closed analytic subset Lk^{Wj) C Wj, and Sk,._-i^iWj) C Wj. By 
forming the disjoint union of those desingularizations and M — M^ one obtains a desin- 
gularization of Mo say 

TTo : M^ ^ M, where 7ro(M^) = Mq. 

Then (|2.1|) restricted to Mo, that is (a, /)|Mo, can be naturally lifted, using the puUback 
operation, to an IDE (ai, /i) = ttq ((a, /)|Mo) on as follows 

a'iiy)y = a{iTo{y))TyTCo{y,y) 

My) = fMy))- 
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We should remark at this point (see also the paragraph Some notation and operations 
with IDE, in section |21) that in the present paper a tangent vector {x, x) to Mq at a point 
X G Mo is a vector {x, x) G T^M such there is an analytic curve z{t) e Mq, say defined 
for t G (—5, S), such that z{0) = x and the tangent vector to z{t) at t = as a curve in 
M coincides with {x,x). In particular, if y{t) is a given analytic curve in then 

TyMy,y) = — -f. — 

is a tangent vector to Mq at ttq (y(0)) . 

Note that is a manifold of dimension dim = dim Mq < ci. 

Complete desingularization of (a, /) in a finite number of steps. Now we 

repeat the process for the IDE (ai, /i) with domain and range F, proceeding as we 
did before with the system (a, /) with domain M and range F. We obtain a decomposition 

= U M| U Ml 

We know that there is no solution to the LAS system 

ai{y)y = fi{y) 

for y G Ml. We also know that there is solution to the same LAS system for y G Mg, 
moreover, (ai,/i)|M2 is an IDE of constant rank. Now we desingularize Mq 

TTl : ^ M\ TTi{M^) = M^ 

and repeat the process. Finally, we obtain a finite sequence of manifolds and maps 

M« "-^^ M"'' ""^'...^M'^M, 

where 7ro(Mi) = Mq, 7ri(M2) = and in general Hi{M'+^) = M*, for i = 0, - 1, 
where we have written M^ = M to unify the notation. 

We have obtained a finite recursive procedure that reduces the problem to a finite number 
of IDE of constant rank, namely, the IDE of constant rank (oj, /j)|M2, for i = 0,1, q, 
where we have written (ao, /o) = (a, /), to unify the notation. We will call this a desin- 
gularization process and the sequence of maps tTj and IDE (a^+i, /j+i), z = 0, g — 1 a 
desingularization of {a, f). 
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The collection (02, /2) of IDE (a^, /fc)|M|, k = 0, g, defines a single IDE of constant 
rank in the disjoint union M2 = Llfc=o "^2 > have said at the beginning of the para- 

graph Description of the algorithm. We have a natural projection 712 : M2 — »■ M. This 
IDE (02, /2) with domain M2 and range F is called the desingularizing IDE. 

Remark. As we have said before the range F remains the same throughout the applica- 
tion of the algorithm. However, in practice it is sometimes convenient to apply theorem 
12.11 which may imply a change of F, to simplify calculations. 

4 The Main Results 

In this section we will show in which precise sense the solutions to the desingularizing 
system {0,2, /2) of a given analytic IDE (a, /) are related to the solutions to (a, /). It is 
clear that in certain simple examples of IDE one can show in a more or less direct way 
that solutions to the desingularizing system project via 712 onto solutions to the given 
IDE, and also that solutions to the IDE are such projections. For instance, if one is 
interested in the local behavior of solutions near a singular point of Mq it is sometimes 
enough to use a simple blow-up to desingularize Mq at that point, and one can show in 
certain cases in a simple and useful way the relationship between solutions to the given 
system and solutions to the system obtained by blow-up, see @]. However, in this paper 
we want a more precise and general global result showing that a curve, belonging to a 
certain convenient class of curves, is a solution to a given IDE if and only if it is essen- 
tially the projection via the map 7^2 of a solution to the desingularizing system belonging 
to the same class of curves. This is important, for instance, if one is interested in global 
aspects of solutions, like extension of solutions, or in a description of the family of all 
solutions. In order to be able to use the theory of subanalytic sets we need to define 
carefully a convenient class of curves, which we will do in the next paragraph. 

From now on, we will often use the theorem 6.1 of [^01 that a subanalytic subset of 
dimension 1 of an analytic manifold M is a semianalytic subset. We will also often use 
the fact that the image of a relatively compact subanalytic subset under a subanalytic 
map is a subanalytic subset, see immediately after definition 3.2. Using this we can 
deduce that the image of a relatively compact subanalytic subset under an analytic map 
is a subanalytic subset, which will be also useful for us. 

In what follows, we will work with several types of intervals, like (tq, Ti), [tq, Ti), (tq, Ti] 
or [To,ri]. We will usually assume that tq and ri are real numbers. However, some 
definitions and results are valid also for the case in which the open end of an interval is 
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±oo, that is, for intervals (—00, Ti); (ro,+oo); (—00, +00); [ro,+oo); (— oo,ri]. 

The notion of an as-curve. Inspired by jHHj, definition 3.2, we will define 

Definition 4.1 A subanalytic curve x : (toi^i) ^ M is a suhanalytic map, that is, a 
map such that graph x C M x M zs a suhanalytic subset. We define the notion of a 
suhanalytic curve x : [to^i) ^ M-, x : (to^i] ^ M or x : [to,^i] ^ M in a similar way. 

In definition 14. II since dim(graphx) = 1 we have that graph x is a semianalytic set. 

Lemma 4.2 (a) Let x : [toj^i) M he a continuous suhanalytic curve whose graph 
is a relatively compact suhset. Then there is a uniquely defined continuous suhanalytic 
extension x : [to;i^i] M. A similar result holds for suhanalytic curves x : (to^ ''^i] ~^ M 
or X : (to, ti) M. 

(b) Let X : [to,ti) M be a continuous subanalytic curve whose graph is not a relatively 
compact suhset. Then graph x is closed. A similar result holds for subanalytic curves 
X : (to, ti] ^ M or X : (to, ti) M. 

Proof. First we shall prove (a). We need to show first that the limit of x{t) a.st ^ t^ ex- 
ists. Using corollary 2.8 of jHOl we can deduce that the closure G oiG = graph x in R x M 
is subanalytic and compact, and then also G fl ({ti} x M) is subanalytic and compact. 
Since Gfl ({ti} x M) is nonempty let xi G Gfl ({ti} x M) be given. One can choose local 
coordinates at xi, and for any small e > the set = Gn((to, ti) x B^{xi)) is relatively 
compact. We can easily deduce from theorem 3.14 of ^Ij that any relatively compact 
subanalytic set has a finite number of connected components. Let Cj, i = I, ...,n(e) be 
the connected components of G^. It is not difficult to see that each connected component 
Ci is of the type Ci = graph (x|(aj, , i = 1, ...,n(e). We can assume without loss of 
generality that Pi < a^+i, i = l,...,n(e) — 1. Since Xi is a limit point of G^ we must 
have j3n(^) = ti, which implies that x(t) G B^{xi) for all t G («n(e), /5n(e)), as we wanted 
to prove. The fact that the continuous extension x : [to, ti] — > M is a subanalytic curve 
follows from corollary 2.8 of [30j. The rest of the proof of (a) can be performed in a sim- 
ilar way. Now we will prove (6). If graph x is not relatively compact then GCi ({ti} x M) 
must be empty, otherwise we can proceed as in the proof of (a) and we can conclude 
that the limit of x(t) as t ^ x^ exists and then, one can show that G is relatively 
compact. Since all the limit points of G not in G must belong to G fl ({ti} x M) we 
have that G is closed. The rest of the proof of (b) can be performed in a similar way. ■ 

In order to define a convenient class of curves to solve a given IDE we introduce the 
following notion 
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Definition 4.3 (a) An analytic- semianaly tic- curve x{t), t G {tQ,ti), in M, where M is 
a given manifold, is an analytic map x : (to,ti) ^ M which is also a subanalytic curve, 
that is, such that graph a; is a semianalytic subset o/R x M. We will often call such an 
analytic- semianalytic- curve in M simply an as-curve in M. 

(6) An analytic-semianalytic-curve (or as-curve) x(t), t G [to,ti), {t G (to,ti],t G [t^^i]) 
in M, where M is a given manifold, is a continuous map x : [to;^i) M, (respectively, 
X : {to,ti] M, X : [to^^i] M), which is also a subanalytic curve, that is, such that 
graphx is a semianalytic subset o/R x M, and x\{to,ti) is an as-curve in M. 

For instance, x = y/t, t G (0,c), (or t G [0, c), t G (0, c], t G [0,c]), with c > 0, are 
as-curves in R. On the other hand, x = tsm{7i/t), t G (0,c), with c > 0, is not an 
as-curve in R, but x = tsm{7c/t), t G {6, c), with < 6 < c, is an as-curve in R. 

Next, we will give some lemmas where some basic properties of as-curves, that we need 
to prove our main results, are proved. 

Lemma 4.4 (a) Let x : (toi^i) ^ N be a given analytic map, where N is a given 
manifold. Then any map x|(to,^i) '■ (^O)^i) ~^ where and ii are such that to < 
to < ii < ti, is an as-curve in N. In a similar way, any map x\[tQ,ii) : [to,ii) — > A^, 
x|(to,ii] : (^O)^i] ^ ^) or x|[to,ti] : [io,ii] N, with and ii as before is an as-curve 
in N. 

(6) Let X : {to,ti) N be an as-curve and assume that graph x is a relatively compact 
subset o/R X A^. Then there is a unique continuous extension x : [to;^i] ^ M which is 
an as-curve. 

Proof. Let us prove (a). We have that graph(x| {to, ti)) is a semianalytic subset of R x 
defined as {(t, x) G (to, ti) x N : to < t < ii , x = x{t)}. The rest (a) can be proved in a 
similar way. To prove (6) we simply apply lemma ■ 

Lemma 4.5 Let x : [to,ti) ^ M be a subanalytic map and assume that x is continuous 
at to. Then there exists t2 G (tc^i) such that x\[tQ,t2] is an as-curve. A similar result 
holds for a subanalytic map x : {to,ti] — M continuous at ti, that is, there exists 
t2 G ito,ti) such that x\[t2,ti] is an as-curve. 

Proof. If X is a constant the result follows immediately. Let us assume that x is not 
a constant. It is easy to see that graph(x| [to, ^2]) is a semianalytic subset of R x M of 
dimension 1, for every ^2 ^ (^O)^i)- We can assume without loss of generality, using, for 
instance, Whitney embedding theorem, that M C ?7 is an analytic submanifold of U, 
where ?7 is a real finite dimensional vector space. Let pi : M. x U ^ R be the projection 
onto the first factor. Continuity of x at to implies that one can choose b G (to,ti), such 
that graph(x| [to, &]) is a relatively compact subanalytic subset of R x f/. According to 
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lemma 3.4 of [301 we have that graph(a;| [to, &]) is a finite union of connected smooth 
semianalytic subsets A such that, for each A, rank(pi|A) is constant. It is not difficult 
to see that each A has dimension or 1 and that rank(pi|y4) is or 1. Moreover, we can 
see that there must be an A, say A = Aq, such that rank(pi|74o) = 1, Pi(^o) = (^0,^2), 
for some t2 G (to,b], and pi{Ao) = [^0,^2]. From this we can easily deduce that x\{tQ,t2) 
is an as-curve and moreover, using lema 14.41 (6), that a;|[to,t2] is an as-curve. The rest 
of the proof can be performed in a similar way. ■ 

Lemma 4.6 (a) Let x : [to,^i] — M, x : [to,^i) M or x : {tQ,ti\ M he an as-curve 
in M. Then x\{tQ^ti) is an as-curve in M. 

(b) Let X : {to,ti) ^ M be an as-curve in M and assume that there is a continuous 
extension x : [to;^i] M, x : [tQ,ti) ^ M or x : (to)^i] ~^ M. Then x is an as-curve in 
M. 

(c) Let f : M ^ N be a given analytic map. Let x(t), t G [to^i) it G (to,ti],t G [to^i]) 
t G (to,ti)) be an as-curve in M. Then f {x(t)), t G [t^^i], is an as-curve in N, for 
each [to,ii] C [to,ti), (respectively, [toji] ^ {to,ti], [ioji] C [to,ti], [toJi] ^ (to,ti)). If 
graph X is relatively compact then / , t G [to, ti), (t G (to, ti], t G [to, ti], t G (to, ti)) 
is an as-curve and graph(x o /) is relatively compact. 

Proof. Part (a) is an immediate consequence of the definitions. Part (6) follows eas- 
ily using corollary 2.8 of jSO]- To prove (c) observe first that fox: (to,ti) — is 
an analytic map. Since graph(x| [to, ti]) is a semianalytic compact subset of M x M of 
dimension 1 we have that graph(/ o x|[to,ti]) = (Ik x /) (graph(a;| [to, ti]) , taking into 
account that 1]rx/:MxM— >-MxA^isan analytic map, is also semianalytic, and a 
similar proof can be given for the case of the intervals (to, ti], [to, ti], (to, ti). If graph x is 
relatively compact then according to lemma we have a continuous extension x which 
is an as-curve and therefore we can apply the first part of (c) to this extension The rest 
of the proof follows easily. ■ 

Lemma 4.7 Let N be a given manifold and let x : [to,ti) ^ N be an as-curve in 
N which is not a constant. Then there exist t2 G (to,ti) such that x((to,t2)) is an 
analytic submanifold which is also a semianalytic subset, x|(to,t2) : (to,t2) x{(tQ,t2)) 
is an analytic diffeomorphism and x|[to, t2] : [to,t2] — > x([to,ti]) is an homeomorphism. 
Similar results hold for an as-curve x : (to, ti] — > A^. 

Proof. Let to < ti < ti, then we have that graph(x| [to, ti]) is a compact semianalytic 
subset of M X of dimension 1. Let Xj, i = 1, n be local analytic coordinates centered 
at x(to) therefore Xi(to) = for z = 1, .., n. Without loss of generality we can assume that 
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Xi(t) is defined for all i = I, --^n and all t G [to;^i]- For some index, say j G {1, 
we must have that Xj{t) is not a constant. We are going to show that there exists 
^2 G such that Xj : [to, ^2] satisfies certain conditions from which the lemma 

follows. Since Xj{t) is the projection on the j-coordinate axis of the curve x\ [to, ii] using 
lemma l^ni (c) we have that graph [to, ti]) C M x M is a semianalytic compact subset 
of dimension 1. The restriction to graph(xj| [to, ti]) of the projection p2 : x ^ 
onto the second factor satisfies p2 it,Xjit)) = for all t G [to,ti]. Then according to 

lemma 3.4 of jSU] graph [to, ti]) is a finite union of connected smooth semianalytic 
subsets A such that rank(p2|^) is constant on A, and it is easy to show that rank(p2|^) 
can take only the values or 1. It is also easy to see that for at least one such A one must 
have that rank(p2|^) = 1, 'P2{.A) = xj ((to, t2)) for some t2 G (to, ti] and P2{A) = p2{A) = 
Xj ([to,t2]) , and therefore that a;j((to, t2)) is an open interval. Moreover Xj|[to,t2] is in- 
jective. and we have that (a;j|(to, t2))~"'^ : Xj ((to,t2)) — *■ (to,t2) is analytic and also its 
graph is a semianalytic subset of M x M and then, because of lemma U3J (6), that there 
is an extension (a;j|[to,t2])~"'^ : a;j([to,t2]) [to^h] which is continuous, and therefore 
a;j|[to,t2] : [to,t2] — > Xj([to,ti]) is an homeomorphism. Let t(s), s G [xj(to),Xj{t2)], be 
the map (a;| [to, t2])7^, in other words, the parameter s represents the coordinate xj. Then 
we have that a;([to, t2]) = {(a^i, : Xi = Xi (t(s)) , s G [xj(to),Xj(ti)] ,i = 1, n} and 

also x((to,t2)) = {{xi,...,Xn) : Xi = Xi(t{s)),s G {xj(to),Xj{ti)) ,i = l,...,n}. /^From 
this we can easily deduce the assertion of the lemma for the case of an as-curve x : 
[to, ti) —>■ N. The case of an as-curve x : (to, ti] N can be proven in an entirely similar 
way. ■ 

Lemma 4.8 Let x{t), t G [to,ti) he an as-curve in N, which is not a constant. Then 
there is a t2 E (to,ti) such that a;([to,t2)) — {^^(to)} is nonempty and locally connected 
at x(to), more precisely, x ((to,t3)) is a neighborhood of xit^) in x ([to,t2)) — {x(to)}, for 
allts G (to,t2). Moreover, x : [to,t2] — >■ a;([to,t2]) is an homeomorphism, x((to,t2)) is an 
analytic submanifold which is a semianalytic subset of N and x : (to,t2) x{{to,t2)) is 
an analytic diffeomorphism. A similar result holds for an as-curve x{t), t G (to,ti], in 
M. 

Proof. We can show using lemma UTI that there exists t2 G (to,ti) such that a;|[to,t2] is 
injective, and that x|[to,t2] is an homeomorphism onto its image and that x : (to,t2) — »• 
X ((to,t2)) is an analytic diffeomorphism where x ((to,t2)) is an analytic submanifold. In 
particular, we have that x{t) 7^ a;(to) for all t G (to, t2]. Let r > small be given. Working 
in local analytic coordinates centered at a;(to) we can show that continuity of x{t) implies 
that there exists t^ G (to,t2] such that x {[to,tr]) C Br (a;(to)) • It can be easily shown 
that x([to,tr)) — {x(to)} = x((to,tr)) is connected. Moreover, for each s G [to,tr) there 
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is an open ball Bs {x{s)) C Br (a;(io)) > with S = S{s), such that Bs (x(s))nx ([tr,i^^2]) = 0- 
Then the open set 

W= [j Bs{x{s)) 

satisfies W Br {x{to)) and W H {x ([^0,^2)) — {^ito)}) = ^ {(to,tr)) ■ This shows that 
X ([^0)^2)) — {^^(^o)} is nonempty and locally connected at a;(to) and also that x ((to^^r)) 
is a neighborhood of x{to) in x ([^0)^2)) — {^^(^o)}- Now for each ^3 G {to,tr] take 

W,,= U B-sixis)) 

where S = Sit^, s) and Bg {x{s)) satisfies B^ {x{s)) C Br {x{0)) and x ([ta, t2])nBg {x{s)) = 
0. Then x ((to, ^3)) = Wt^ fl (x ([to, ^2)) — {x{to)}) . This shows that x {{to, t^)) is a neigh- 
borhood of x(to) in X ([to, ^2)) — {x{to)}- Then the first case of the lemma is proved. The 
case of an as-curve x{t), t G (to,ti], in M can be proven in an entirely similar way. ■ 

Inspired by lemma 6.3 of [20] we will prove the following result about the image and 
reparametrization of an as-curve. 

Lemma 4.9 Let x{t) G M, t G [to,ti), be an as-curve in M. Then there is an as- 
curve z{s), s G (so — 61, So + ^2), in M, for some 61,62 > such that z{so) = x{to), 
z{[so. So + 62)) = x{[to, ^2)) for some t2 G (0, ti) and we also have that t2 and sq + 62 are 
such that x\[to,t2) and z\[so, so+62) are homeomorphisms onto z{[so, 80-1-62)) = a;([to, ^2)), 
a;|(to,t2) and z\{so,so + 62) are analytic difeomorphisms onto z{{so,so-\-62)) = x{{to,t2)) 
which is an analytic submanifold which is also a semianalytic subset and, moreover, 
(x|(to, ^2))""^ o {z\{so,So + ^2)) : (-50, -50 + 62) ^ M is an as-curve in R which is an 
analytic diffeomorfism onto its image (to, ^2)- Moreover, t2 can be chosen such that for 
each ts G (to,t2], a;((to,t3)) is a neighborhood of x{to) in x([to,t2)) — {x(to)}. Similar 
results hold for as-curves x{t) in M, where t G (to,ti]. 

Proof. Using lemma and also lemma 6.3 of [201, "^^^ conclude that there is an 
as-curve z{s), s G (sq — 61, so + ^i), in M, for some 61,62 > such that z{so) = x{to), 
z{[so,so + 62)) = x([to,t2)) for some t2 G (0,ti) satisfying all the conditions stated in 
lemma Since we can also apply lemma 1^31 and lemma to z{s) we can also deduce 
that t2 and 62 can be chosen such that x|[to,t2) and 2;|[so,So + ^2) are homeomorphisms 
onto z{[so, 80-1-62)) = x([to,t2)), x|(to,t2) and z| (50,^0 + ^2) are analytic diffeomorphisms 
onto z{{so. So + 52)) = a^((to,t2)) which is an analytic submanifold which is also a semi- 
analytic subset and, moreover, (2;|(to, ta))"^ o (x|(to,t2)) : (to,t2) — > M is an as-curve in 
M which is an analytic diffeomorfism onto its image (so, so + 52)- The rest of the proof 
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can be performed in a similar way. 



The notion of a pas-curve and of a Ics-curve. By definition, a continuous curve 
in a manifold M, x{t), t G [a, h) such that there is a partition a = to <,...,< tm = b such 
that the restrictions x|[tj,tj+i], i = 0, ....,m — 2, x\[tm-i,tm) are as-curves in M is called 
a continuous piecewise-as-curve in M (or continuous pas-curve in M). Each restriction 
x\[ti,ti+i), x\{ti,ti+i], x\{ti,ti+i), x\[ti,ti+i], i = 0, ....,m-2, x\[trn-i,b), x\{tm-i,b) is 
called an as-piece of x. We define the notion of a continuous piecewise-as-curve (or con- 
tinuous pas-curve) x{t), t G (a, 6], t G [a, 6] or t G (a, 6), in M, in a similar way. 

We shall also introduce the notion of a (not necessarily continuous) piecewise-as-curve 
in M (or pas-curve in M) by eliminating from the previous definition the condition of 
continuity at the points x(ti), i = I, ...,m — 1, belonging to two consecutive as-pieces, 
and replacing it by the weaker condition of left or right continuity. More precisely, a 
pas-curve x : [a,b) —>■ M in M is defined by the condition that there is a partition 
a = to <,...,< tm = b such that, for each i = 0, ...,m — 2 such that tj ^ tj+i, the 
restriction x\{ti,ti+i) is an as-curve in M having compact graph and x\{tm_i,b) is an 
as-curve in M having a continuous extension x\[t^_-^,b) obtained by taking the limit at 
tm-i on the right. Besides, we require that x be left continuous or right continuous at 
each ti, i = l,...,m — 1. We define the notion of a piecewise-as-curve (or pas-curve) x(t), 
t & {a,b], t & [a, 6] or t G (a, b), in M, in a similar way. 

We see immediately from this definition that any as-curve is a pas-curve and, moreover, 
that any continuous pas-curve is a pas-curve, in a natural way. Using lemma we can 
conclude that, given a pas-curve x : [a,b) M as before, for each i = 0, ...,m — 2 there 
is a well defined as-curve denoted a;| [t^, tj^^^], which, for each i = 0, ...,m — 2 such that 
ti 7^ tj+i, is the uniquely defined continuous extension of the restriction x|(tj,tj+i). We 
can easily see from lemma that if the graph of x\[t^_i,b) is not relatively compact 
then it must be closed. Similar statements hold for pas-curves x{t), t E (a, 6], t G [a,b] 
or t G (a, 6), in M. 

Given a pas-curve x{t), t G [a,b), as above, we have that some of the restrictions 
x\[ti,ti+i), x|(tj,tj+i], x|(tj,tj+i), for each i = 0,...,m — 2, and some of the 

restrictions x\[tm-i,b) \(tm-i,b), satisfy the required continuity condition at the closed 
end of the interval and then they are as-curves. Each one of those restrictions which is an 
as-curve is called an as-piece of x. If the graph of x is relatively compact we will call each 
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as-curve a;| [t^, t"^]^], i = 0, ...,m — 2, x\ [t^_i, b~] an as-piece closure of x. If graph x is not 
relatively compact then the as-piece closures are x\ [tf, t~^_-^], i = 0, ...,m — 2, x\ [t^_i, h). 
We must observe that an as-piece closure is an as-curve but it is not always an as-piece 
of X. By definition, the open as-pieces are x\(ti, U+i), i = 0, ...,m — 2, x\{tm-i,b). Similar 
definitions hold for curves x{t), where t & {a,b], t E [a,b], t E (a, b). 

We need the following definition. A pas-curve ?/ in M is a refinement of a pas-curve z in 
M if graph?/ = graphs and the graph of each as-piece of y is contained in the graph of 
some as-piece of z. 

We introduce the notation C"'^([a, b), M) to denote the set of all as-curves x : [a, b) M, 
where M is a given manifold. In a similar way, we define C"'*((a, 6], M), C"^'*((a, 6), M), 
C"^([a,5],M), and also ^^"^([a, 6), M), C''"^((a, 6], M), CP'^^((a, 6), M), C^"^^([a, 6], M). 

Lemma 4.10 (a) Every pas-curve in the manifold M, say x{t), t G [a, 6), (t G (a, 6], 
t G [a,b], t G (a, b) ), is a suhanalytic curve in M. 

(6) Let x{t), t G [a,b), be a suhanalytic curve in M whose graph is a relatively compact 
subset o/ M X M and which is left or right continuous at each t G [a,b). Then there 
is a, not necessarily unique, pas-decomposition of x{t), that is, a partition of [a,b], say 
a = to < ti, < tr = b, such that x\[tf ,t~^i], i = 0, r — 2 and x\ [t^_i, b) are as-curves. 
A similar result holds for suhanalytic curves x{t), t G (a, 6], t G [a, b] or t E (a, b). 

Proof. The proof of (a) is a consequence of the fact that each as-piece of a pas-curve 
has a semianalytic graph and also the fact that a finite union of semianalytic subsets is 
a semianalytic subset. To prove (6) we observe first that, because of lemma there 
is a continuous extension, x{t), t G [a,b], of x{t). Also lemma implies that for each 
t G [a,b] the right and left limits x{t~), x{t~^) exists. Then, by lemma HTSl for each 
t G (a, 6) there exists e > such that x\[t~^, {t -\- e)~], x\[{t + e)~^,t~] are as-curves, and 
also x\ [a+, (a + e)"], x\[{b — e)~^, b~] are as-curves for some e > 0. Then the proof follows 
by a standard compactness argument. ■ 

Since we want to work within the subanalytic category, it seems that a good choice for a 
class of possible solutions to a given IDE would be the class of subanalytic curves in M 
whose graph is a relatively compact subset and which are left or right continuous at each 
point. The decomposition, predicted in lemma E.lUl of such a curve in as-pieces, which 
gives its structure as a pas-curve in M, which is not unique, is what will enable us to 
define, in the next paragraph, how such curves can be interpreted as being pas-solutions 
to a given IDE. In other words, one convenient class of possible solutions to a given IDE 
would be the class of subanalytic curves having a relatively compact graph, interpreted 
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as pas-curves, according to lemma E.lUI As another choice for a class of possible solutions 
one could choose the class of all pas-curves in M. We are going to work also with a bigger 
class of curves, namely the class of curves that are locally-compact-subanalytic, or Ics- 
curves, in the sense that their restriction to any compact subinterval is a subanalytic 
curve whose graph is compact and which are left or right continuous at each point. More 
precisely, we define 

Definition 4.11 Let M be a given manifold. Then we define the following classes of 
locally-compact-subanalytic curves, also called Ics-curves: 

(a) C'"'^^ {[tQ,ti), M) , where ti G M orti = +00 is the set of all functions x : [to;^i) M 
such that x\[to,t2] is a pas-curve, for each t2 G (to,^i)- 

(b) C^"' {{to, ti], M) , where to e K or to = —00 is the set of all functions x : (to, ti] — > M 
such that x|[t2,ti] is a pas-curve, for each t2 G (to,ti)- 

(c) C''^* ((to, ti), M) , where to & M. or to = —00 and ti G M or ti = +00 is the set of all 
functions x : (to,ti) M such that x|[to,t2] is a pas-curve, for each [to,ti] C (to,ti). 
{d) C''^^ {[to, ti], M) , where to G M and ti is the set of all functions x : [to, ti] — > M 
such that X is a pas- curve. 

We may define the notion of an as-piece and also the notion of a refinement of a given 
Ics-curve in a similar way as we did in the case of a pas-curve. 

Remark, (i) With the notation of definition 14.111 'we have, respectively for (a), (6), (c) 
and {d), that each one of the graphs, graph(a;| [to, t2]), graph(a;| [t2, ti]), graph(x| [to, ti]) 
and graph(x| [to, ti]) is compact. 

(ii) In view of lemma E. 101 we may replace pas-curve by subanalytic curve which is left or 
right continuous at each point in definit ion 14 . 1 1 1 and we will obtain an equivalent defini- 
tion. The decomposition of each subanalytic curve a;| [to, t2], x\ [t2, ti], x| [to, ti] or a:| [to, ti] 
as a pas-curve in M is not unique. For a given Ics-curve in M there is a decomposition 
in at most a countable number of as-pieces, which is not unique. 

Solutions to IDE: as-solutions, pas-solutions and Ics-solutions. Now we intro- 
duce the notion of solution to a given IDE which is convenient for the purposes of the 
present paper. 

Definition 4.12 (a) An as-solution x{t), t G [to,ti) in M to a given IDE {a, f) in M is 

an as-curve in M which satisfies {a, f) for all t G (to, ti), that is, a {x{t)) x{t) = f {x{t)) , 
for all t G (to,ti). Similar statements hold for as-solutions x{t) in M where t G (to,ti], 
t G (to, ti) or t G [to, ti] with to 7^ ti. 

(6) A pas-solution in M to a given IDE is a pas-curve in M such that each open nonempty 
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as-piece is an as-solution in M. 

(c) A Ics-solution in M to a given IDE is a curve in M belonging to C^^" ([to,ti),M) , 
C""' ((to, ti], M) , C^^' ((to, ti), M) or C^^' ([to, ti], M) swc/i that each restriction to a com- 
pact suhinterval has a pas- decomposition which is a pas-solution in M. 

Lifted and projected solutions. We have the following result. 

Lemma 4.13 (a) Let y{t), t G [to,ti), (t G (to,^i], ) be a given as-solution in to the 
system {ak, fk) described in the previous section, for some k = 1,2, ...g. Then for each 
t2 G (to,^i), y(t) is projected via Hk-i into an as-solution x(t) to the system (a^-i, fk-i), 
x{t) = TTfe^i {y{t)) , t G [to, ^2], (respectively, t G [t2,ti] ), in M^'^ . 

ip) Assume that we have an as-solution y{t), t G [to,ti), (t G (to,ti] ), in , for some 
k = 1,2, ...g, to the system {ak, fk) described in the previous section. Then for each 
s = 0, k — 1, y(t) is projected via Tig o ... o Tr^-i into an as-solution x{t) to the system 
{as,fs), x{t) = TT, o ... o Tik-i (y(t)) , t G [to,t2], (respectively, t G [t2,ti], ), in M", for 
each t2 G (to, ti). 

Proof. Part (a) is easy to prove using the fact that tt^-i is an analytic map, and also 
lemma (c). Part (6) follows using (a). ■ 

^From the previous lemma we can deduce that if y(t), t G [to,ti) is an as-solution to 
[ttk, fk) in M'', for some k = 1, ...,q and x(t) = vr^o, ...,7ik {y(t)) is an as-curve in Af, for 
some s = l,...,k, then x(t), t G [to,ti) is an as-solution to (as,/^). We will call y{t) a 
lifted as-solution of x(t), and x{t) the projected as-solution of y{t). A similar definition 
holds for the case of as-solutions ?/(t), t G (to,ti], t G (to,ti), t G [to,ti]. The notions 
of a projected and lifted pas-solution are defined as follows. A pas-solution to {ak,fk), 
say y{t) G M'', t G [to,ti) is a lifted pas-solution of a pas-solution x(t) to {ttg, fs) if 
each as-piece of y(t) is projected, via the composition tt., o ... o tt^^i, onto an as-piece 
of x(t) = vTs o ... o 7[k-iy{t), which is then an as-solution to {as,fs)- We call x(t) the 
projected pas-solution of y(t). A similar definition holds for the case of pas-solutions 
y{t), t G (to,ti], t G (to,ti), t G [to,ti]. 

Main results. 

Theorem 4.14 (a) Let y(t), t G [to,ti) (respectively, t G (to,ti],j be an as-solution to 
{ak,fk) in M'', k = 1, ...,q. Then x(t) = iTk-i {y{t)) , t e [to,t2] (respectively, t e [t2,ti]) 
is an as-solution to {ak-i, fk-i) in M^^^, for each t2 G (to,ti). The previous statement 
holds true if we replace an as-solution by a pas-solution or a Ics-solution. 
(b) Ifx(t), t G [to,ti) (respectively, t G {to,ti]) is an as-solution to (ctfc-i, /fe-i) in M''~^ 
such that x{t) G Mq~^, t G [to,ti) (respectively, t G (to,ti]j, k = 1, then there exists 
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h ^ (i^O)^i) fl^'^ lifted as-solution y(t), t G [to)'^2] (respectively, t G [^21^1]^' o/x|[to,i2] 
(respectively, x\[t2,t^) to {ak, fk) in , in particular, x{t) = Tik-i {y(t)) , t G [^07^2] 
(respectively, t G [^2,^1]-^' 

(c) Let x{t), t G [to,ti] be an as-solution to {ak-i, fk-i) in M^~^ such that x{t) G Mq~^, 
t G [toj^i]; k = l,...,q. Then there is a lifted pas-solution of a refinement of x, say 
y(t), t G [to)^i] ^0 {ak, fk) in M^, in particular, x(t) = ^k-i {y(t)) , t G [^o^^i]- More- 
over, if C 7r^_^-,^(x([to, ^i])), k = l,...,g, is a compact semianalytic subset such that 
TTk-i{C) = x{[to,ti]) then there is a lifted pas-solution of a refinement ofx, say y{t) G C, 
t G [to,ti], to {ak, fk) in , in particular, x{t) = TCk^i {y{t)) , t G [to,^i]- 
{d) Let x{t), t G [to,ti] be a pas-solution to (flfc-i, fk~i) in M^~^ such that x{t) G M^^^ , 
t G [to7^i]) k = l,...,q. Then there is a lifted pas-solution of a refinement of x, say 
y{t), t G [to,ti\ to {ttk, fk) in M'^, in particular, x(t) = TTk-i {y(t)) , t G [to;''^i]- More- 
over, if C C 7i'^}^-^^(x{[to,ti\)), k = l,...,q, is a compact semianalytic subset such that 
7rfc_i(C) = x{[to,ti]), then there is a lifted pas- solution of a refinement ofx, say y{t) G C, 
t G [to,ii], to [ak, fk) in , in particular, x{t) = Tr^-i {y{t)) , t G [to,ti]. 
(e) Let x{t), t G [to, ^i), (t G {Iq, ti], t G [to, ti], t G (to, ti) ) be a Ics-solution to (flfc-i, fk-i) 
in M^~^, such that x{t) G Mq~^, t G [to,ti), (respectively, t G (to,ti], t G [to,ti], 
t G (to,ti)j, k = l,...,q. Then there is a lifted Ics-solution of a refinement of x, say 
y{t), t G [to,ti), (respectively t G (to,ti], t G [to,ti] or t e (to,ti)j to {ak, fk) in , in 
particular, x(t) = nk-i {y(t)) , t G [to,ti), (t G (to,ti], t G [to,ti], t G (to,ti)j. More- 
over, ifCC 7r^_^^(x([to,ti))), (respectively, C C TT^\{x{{to,ti])), C C 7r^i^(x([to, ti])), 
C C 7r^^-^(x((to, ti)))^, k = l,...,q, is a subset such that vrfc_i(C) = x([to,ti)), (respec- 
tively, nk-i{C) = x((to,ti]), Tik-iiC) = x([to,ti]), 7ik-i{C) = x{{to,ti))), and, besides, 
C n 7r^i;^(a;([to, ti])) is compact subanalytic for each compact subinterval [to,ti] C [to,ti) 
(respectively, [to,ti] C (to,ti], [to,ti] C [to,ti], [to,ti] C (to,ti)j, then there is a lifted 
Ics-solution of a refinement of x, say y(t) E C, t E [to,ti), (respectively t G (to,ti], 
t G [to,ti], t G (to,ti)) to {ak,fk) in M^, in particular, x(t) = vr^.i (?/(t)) , t G [to,ti), 
(t G (to,ti], t G [to,ti], t G (to,ti);. 

Proof. We are going to give a detailed proof of the case k = 1 only, since the cases 
k = 2,...,q can be proven in an entirely similar way. Part (a) is an easy consequence 
of lemma 14.131 In order to prove (6) , (c) , {d) and (e) we are going to prove first sev- 
eral facts, namely, {i), {ii), (Hi), (iv) and (v) below. These facts will be proven under 
the assumption that x(t) G Mo, t G [to,ti], is an as-solution to (a,/) in M, the curve 
x{t) is simple, that is, x{a) 7^ x{b) for all a,h E [to,ti] such that a ^ h, and moreover, 
X : [to,ti] — s> x([to,ti]) is an homeomorphism and also x((to,ti)) is an analytic subman- 
ifold and x|(to,ti) : (to,ti) — >■ x((to,ti)) is an analytic diffeomorphism. We can assume 
without loss of generality (for instance, using Whitney embedding theorem) that M C f/ 
and C V are analytic submanifolds of U and V, where U and V are real finite di- 
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mensional vector spaces. 



{i) The map ttq can be described as the restriction p| (graph ttq) to graph ttq (1 V x U 
of the projection onto the second factor p : V x U U. Since p| (graph ttq) is a proper 
analytic map we have that (p| (graph 7ro))~^ {x ([^0,^1])) is a compact semianalytic subset 
of V X U, therefore, according to lemma 3.4 of [301 it is a finite union of relatively 
compact connected smooth semianalytic subsets A such that for each A rank(p|y4) is 
constant on A. Since dim (a;([to, ^i]) = 1, it is easy to see that for each A Ta.nk{p\A) 
is or 1, and moreover, x{[to,ti]) is the union of those p{A) such that p\A has rank 
1 and therefore p{A) is not a point. We observe that if there is a compact semi- 
analiytic subset C C n^-^ {x{[to,ti])) such that vro(C) = x([to,ii]) then we have that 
Ci = (p| graph ttq)^^ {x{[tQ,ti])) H {C x U) is compact and semianalytic. According to 
lemma 3.4 of [201, Ci is a finite union of relatively compact connected smooth semian- 
alytic subsets A having the same properties as before. 

[a) Let a G [to)^i] be fixed. Then there is an A, say A = Aq, such that x{a) G p{Aq). 
Since Aq is connected and compact and the curve x{t), t G [to, ti], is simple, we have that 
p{Aq) is homeomorphic to a closed interval (possibly of zero length) via the curve x, say 
p{Ao) = a; ([01,02]) , where [01,02] C [to,ti]. We can assume without loss of generality 
that p{Ao) is homeomorphic to a closed interval of nonzero length. 

[in) By using theorem 6.10 of (20] we can see that for given points G p~^ {x{ai)) hAq, 
i = 1, 2, so in particular, we have qi ^ q2, there is a continuous semianalytic curve w{s), 
s G [so,si], in Aq such that w{si) = qi+i, ^ = 0,1. Then w([so,si]) is a compact semi- 
analytic subset of V X U of dimension 1 and we have p {w{[sq, si])) = x{[ai, 02])]. Using 
lemma 3.4 of [20] we see that since w([so,si]) has dimension 1 it is a finite union of 
relatively compact connected smooth semianalytic subsets B of dimension less or equal 
than 1, such that the restriction of the projection p\B has constant rank of value or 1. 
Since each B is connected and compact we have that p (i?) is homeomorphic to a closed 
interval, say p [B) = o: ([o^, 6^]) • For at least some B such that p [B) is not a point one 
must have that p (S) contains the point x{a). Observe that p{B) = x {{aB,bB)) is an 
analytic submanifold which is a subanalytic subset and that p\B : B x((ob,6b)) is 
an analytic diffeomorphism. This gives, in particular, a parametrization of the analytic 
submanifold B with the parameter t, namely, z{t) = {p\B) ^ {x{t)) , t G {aB,bB)- We 
can give a definition of graph (^{p\B)~^ ox) as a subanalytic subset of M x x M as 
follows. We have graph [{p\B)~^ o x) = {{t, z) e R x x M : z e B,p{z) = x, {t, x) G 
graph(x|(oB, &b))}, which defines graph (^{p\B)^^ o x) by subanalytic conditions, since B 
is a semianalytic subset of x M, graph x is a subanalytic subset of R x M, by def- 
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inition, and graph(x|(aB, 6^)) = E M. x M : (t,x) E graphx, t G (aB,&B)}- Then 

using lemma 14.21 we can deduce that there is a uniquely defined continuous extension, 
which we will call z hy a. slight abuse of notation, t G [05,6^], whose image is 
the semianalytic subset B = z{[aB,bB]), which is an as-curve in x M. We have, in 
particular, that z(t), t G [aB,^^] is an as-curve such that p{z{t)) = x{t), t G [0^,65]. It 
is clear that the extension z(t), t G [05,^5] is given hy {p\B) ox 

{iv) Assume that x{o.) G p{B). We have the as-curve z(t) = {y{t),x{t)) , t G [05,65] 
in graph(7ro), therefore p{z{t)) = x(t), for t G [aBybs], then the curve y(t) satis- 
fies 7io{yit)) = x{t), t G [05,65]. Using this it becomes clear from the definition of 
(ai,/i) that y{t) satisfies the system (ai,/i), for t G {aB,bB)- We have that, for any 
[a — ei, + 62] C [05, 65], where ei, £2 > 0, ?/(t), t G [« — ei, + 62] is an as-curve in M^. In 
fact, this is a direct consequence of lemma (c), since the projection q : x M — > 
is an analytic map. It is clear that y(t) satisfies the system (ai, /i), for t G (a — ei, a + e2). 

(f) Assume now that x{a) G p — p , then a = 05 or a = 65. If a = 05 (re- 
spectively a = 65) we can proceed in a similar way as we did in (iv) and we have an 
as-curve z(t) = {y(t),x(t)) , t E [a,a + €2] (respectively, t E [a — ei,a]) in graph(7ro), 
then, in particular, p{z(t)) = x(t), t E [a, + 62], (respectively, t E [a — ei,a)]. Then the 
curve y{t) = q {z{t)) satisfies ttq {y(t)) = x{t), t E [a,a + €2] (respectively, t E [a — ei,a]) 
and is an as-curve. It is clear that y{t) satisfies the system (ai,/i), for t E {a,a + €2) 
(respectively t E (a — ei, a)). 

We are going to prove (b). If x{t) = x{to) is a constant then it can be lifted to a constant 
curve y{t) = y{to), where y(to) E tiq^ (a;(to)) , which solves the problem, so let us assume 
that x{t) is not a constant. By conveniently lowering the value of ti we can assume 
without loss of generality that x(t) E Mq, t E [to,ti], is an as-solution to (a,/) in M. 
Moreover, by lowering the value of ti if necessary and using lemma ITTI we can assume 
that the curve x(t) is simple, that is x(a) 7^ x{b) for all a,b E [to, ti] such that a b, and 
moreover, that x : [to,ti] — >• x([to,^i]) is an homeomorphism onto x{[to,ti]), x ((to,^i)) is 
an analytic submanifold of M and a;|(to,^i) '■ {to,ti) —>■ x{{to,ti)) is an analytic diffeo- 
morphism. By using (v) with a = to, we must have 05 = to and then the proof of (6) 
follows from (f) by taking a + €2 = ^2- The case of an interval (to;^i] can be proved in 
an entirely similar way. 

We are going to prove (c). As in the proof of (6), the case in which x is constant is 
trivial, so we shall assume that x is not a constant. Using (6) we can conclude that 
for each i E (tQ,ti) there are intervals [i — ei,t\ C [to,ti] and + €2] C [to,ti], where 
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ei,e2 > 0, and lifted as-solutions y^(t), t G [t — ei,t|, and y^(t), t G + €2], in M^, 
of — ei,t\ and t + €2] to (ai,/i). Also if t = to = ti) there is a lifted as- 
solution y^(t) of a;| [to, + ^2] (respectively, of x\[ti — ei,ti],) to (ai,/i). Using this 

and a compactness argument we can conclude that there is a partition of [to^^i], say, 
to = to < ti,...,< ir = ti, such that there are lifted as-solutions yi{t), t G of 
x\[ii, t(+i], to (ai, /i), z = 0, r — 1. Then a lifted pas-solution can be obtained by gluing 
together the lifted solutions yi{t), i = 0, r — 1 into a single pas-solution y(t),t G [to, ti], 
which is always possible, and in more than one way if r > 0. If there is a set C satisfying 
the conditions stated in (c) then we can prove that there is a lifted pas-solution satisfying 
the required conditions proceeding as above and taking into account the comments made 
at the end of (i). 

The proof of [d) can be obtained as a consequence of (c) by conveniently gluing together 
a finite number of lifted pas-solutions each one corresponding to each as-piece of the 
given pas-solution. 

The proof of (e) can be obtained as a consequence of (d) by conveniently gluing together 
conveniently chosen at most countable lifted pas-solutions. ■ 

The next theorem adds some more information and completes the picture of the rela- 
tionship between solutions and lifted solutions to a given IDE. 

Let x{t), t G [a, P] be an as-solution to (a, /) in M*^, for some k = 0, ...,q — 1. If for some 
So G [a, /?] we have x{so) G M| then we must have that there is at most a finite number 
of t G say a < ti < .... < U < P, such that x{ti) G Mq. This is because Mq is 

defined by analytic equations. In this case we will call x a normal as-solution to (a^, fk) 
in M^. If x{t), t G [0.1(3] is a pas-solution to {ak, fk) in M^, for some k = 0,...,q — 1 
such that all its as-pieces are normal we will call x a normal pas-solution to {a^, fk) in M^. 

Let us introduce the following notation. M = \_\'1=qM^, where \_\ means disjoint union, 
vf : M — i> M, where vf (x) = TTk{x) if x G M*^. Now assume that we have a finite collection 
of normal pas-solutions of the following type. For each i = I, ...,r let yi{t), t G [ti,o,^i,i] 
be a normal pas-solution to some of the IDE {ak,fk), k = 0,...,q, say {aki,fki), where 
ki G {0,...,q}, i = l,...,r, so that yi{t) G Mq* only for at most a finite number of 
t G [ti,o,tj,i]. Assume, without loss of generality, that tj^i = ti+i,o, for i = l,...,r — 1, 
so we have a partition ti^o < ti^i = t2,o < ••• < ^r,i of the interval [ti,o,tr,i]- Then the 
collection of curves yi, i = 1, ..,r may be thought of as a single pas-curve in M, say y(t), 
t ^ [^1,05 ^r,i]- We will say that y{t), t G [ti^o, ^r,i] is a normal pas-solution to (a, /) in M. It 
follows from theorem l4.141 (a) that x(t), t G [ti,o, tr,!], given hy x{t) = ttqo ... onk^-iyi(t), 
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t G [tifi, ti^i], i = I, ..,'r, or, equivalently, x(t) = n {y{t)) , t E [ti.o, i^r,i], is a pas-solution 
to (a, /). In this situation we will say that the solution y in M is a normal pas-lifting of 
the solution x in M, or, equivalently, that y is a lifted pas-solution of x, and also that x 
is the pas-projection of y. 

We can use the previous remarks in combination with theorem 14.141 and we obtain the 
following theorem 

Theorem 4.15 Let x{t), t G [«,/?] he an as-solution to {a,f). Then there is a lifted 
normal pas-solution y{t), t G [a,l3] to (a,/) in M of a refinement of x{t). 

Equivalence between an IDE (a, /) in M and (a2, f2) in ^2- Using the notation 
introduced in section |21 we have obviously that M2 is an open subset of M, 112 = vr|M2, 
(a2,/2) = (5,/)|M2. 

Let y{t), t G be a normal as-solution to {ak, fk) in M'^ and let a < ti < .... < 

tr < P, be the set of all t G [a,/3] such that x(ti) G Mq, as explained before the the- 
orem. Then we obtain a well defined collection of as-solutions to {ak, fk)\M2 in M|, 
namely, tj+i), i = l,...,r — 1, y\{a,ti), if a < ti, y\(tr,l3), if tr < p. We will call 
this collection of solutions the collection of solutions in , or, also, in M2, induced 
by the as-solution y. It is clear that by continuous extension, using lemma 14. 2t of each 
y\(ti, tj+i), i = I, r — 1, y\{a, ti), if a < ti, y\{tr, (3), if tj. < /3, we obtain a refinement 
of X, with as-piece closures y\[tf,t^^i], i = l,...,r — 1, if a < ti, y\[t^,P~], if 

tr < p. Of course, by continuity of y we have y\[tt ,t~'^-^] = ?/| [tj, tj+i], i = l,---,r — 1, 
y\[a+,t^] =y\[a,ti], if a < ti, y\[tt,p-] = y\[tr,P], if tr < p. 

Given any normal pas-solution y to {ak, fk) in M^, we shall define the collection of as- 
solutions in M2, or, also, in M2, induced by the normal pas-solution y as being the union 
of the collections of as-solutions induced by all the as-pieces of y. Finally, given any nor- 
mal pas-solution y to (a,/), in M. we shall define the collection of as-solutions in M2 
induced by the normal pas-solution y as being the union of the collections of as-solutions 
induced by all the as-pieces of which y is ultimately composed. 

It is clear that a given normal pas-solution y{t), t G [a, P] to (a, /) in M, say y{t) G M^, 
for some k = 0, g as before, cannot be completely recovered from its induced collection 
of as-solutions to (02, /2) in M2. This is because even if it is known that the value of y(t) 
at some t G [a, /3] is the left or the right limit at t of some of the as-pieces of the induced 
collection of as-curves by y, it is not known which one of those limits, unless one has 
some additional information, like continuity of y at t. Modulo this loss of information. 
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we can consider the IDE (02, /2) with domain M2 and range F and (a, /) with domain M 
and range F as being equivalent. In this sense, the meaning of theorems 14 . 141 and 14 . 1 51 
is that solving a given analytic IDE (a, /) with domain M and range F can be reduced 
essentially to solving a finite collection of IDE of constant rank, namely, (aj,/i)|M2, 
i = 0, g — 1, or, equivalently, the system of constant rank (02, /2) with domain M2 and 
range F. As we have said before it is clear that each (oj, fi)\Ml, z = 0, q — 1 defines an 
analytic family of vector fields and hence an analytic affine distribution of constant rank 
on the manifold M2. This immediately implies that (02, f2) defines an analytic family of 
vector fields and hence an analytic affine distribution of constant rank on the manifold 
M2. Therefore, it also gives rise to an analytic control system, |31j. The usefulness 
of this kind of result relies obviously on the interest that such control systems have in 
several fields. For instance, the well known theorem of Sussmann, jl21; tells us how to 
deal with problems of reachability and observability, once one has a control system. See 
also where this kind of questions are solved in the context of subanalytic sets. 

Reparametrization and extension of solutions. First we shall define the notion 
of a reparametrization in the context of as-curves. A reparametrization of an as-curve 
x(t),t G [tQ,ti), in M, is a change of variables t = r(s), s G [so,si), (s G (si,so]) such 
that r : [sq, Si) — ^ M (respectively, r : (si,So] — > M) is an as-curve in M which is also 
an homeomorphism onto [to^^i) such that r(so) = to and t|(so,Si) : (sq, Si) (tQ,ti) 
(respectively, r : (si,so) (to,ti)) is an analytic diffomorphism. It is easy to prove 
that in this case the composition (x o r)(s), s G [sq, Si) (respectively, s G (si, Sq]) is an 
as-curve in M. Similar definitions and results hold for the case of as-curves x{t) whose 
domain is an interval of the type ito,ti], (to,ti), [^oj^i]- 

Let (a, /) be a given IDE with domain M and range F. By definition, (a, /) is homoge- 
neous if / = 0. It is clear that if x{t), t G [to, ti) is a given as-solution to an homogeneous 
system (a, 0) and t = r(s), s G [sq, Si) (s G (si,so]) is a given reparametrization then 
the curve y{s) = x o r(s), s G [sq, Si) (respectively, s G (si,So]) is also an as-solution 
to (a, 0). More generally, if (a,/) is not necessarily homogeneous then y{s) satisfies 
a {y{s)) y{s) = {dr / ds)f {y{s)) , s G (so,Si), or, using a different and also standard 
notation, a{y{s))y{s) = t{s)f{y{s)), s G (so,si). Similar results hold for the case of 
as-curves whose domain is an interval of the type (to,ti], (to,ti), or [to,ti]. 

Theorem 4.16 (a) Let (a, 0) be a given homogeneous IDE with domain M and range F. 
Let x+{t), t G [to, ti), be an as-solution to (a, 0) in M, which is not a constant. Then there 
is an as-solution z{s), s G (sq—^i, 59+52), in M, for some 61, 62 > 0, satisfying all the con- 
ditions stated in lemma \^~^ More precisely, z{so) = a;+(to), z{[so, 50+^2)) = a;+([to, t2+)) 
for some t2+ G (0,ti) and we also have that t2+ and Sq -|- 62 are such that a:+|[to,t2+) 
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and z\[so, Sq + S2) are homeomorphisms onto z{[so, Sq + S2)) = x+([to, ^2+))) a;+|(to7^2+) 
and z\{so, So + 62) are analytic diffeomorphisms onto z{{so, So + 62)) = x+((to,i2+)) which 
is an analytic submanifold which is also a semianalytic subset and, moreover, the curve 
in M, say t = t{s), given by (a;+|(to, ^2+))"^ ° {z\{so,sq + 62)) : [sq, so + 62) ^ M. is an 
as-curve in M which is an analytic diffeomorfism onto its image (^0,^2+). Moreover, t2+ 
can be chosen such that for each ^3+ G (^0)^2+], x+ ((^0!^3+)) a neighborhood of x+ito) 
in x+ ([to)^2+)) — {^+(to)}- Similar results hold for as-solutions x„(t), t G (to)^i]- 
(6) Let {a, f) be a given IDE with domain M and range F. Let x+{t), t G [to,ti), be 
an as-solution to {a, f) in M which is not a constant. Then there is an as-curve z{s), 
s G (so — Si, So + ^2), in M, for some 61, 62 > 0, not necessarily a solution, satisfying all 
the conditions stated in lemma \J~^ as we have explained in (a), and, besides, the as-curve 
z{s), s G [so, , So + ^2) satisfies the equation a {z{s)) z{s) = i{s)f {z{s)) , s G (sq, Sq + ^2), 
with t{s) > 0, s G (so, si). Similar results hold for as-solutions of the type X-(t), 
t G ito,ti]. 

(c) Let {a, f) be a given IDE with domain M and range F. Let x+{t), t G [to,ti), be an 
as-solution to {a, f) in M which is not a constant and let z{s) be as in (6). Then by 
conveniently diminishing the value of 61 and the value of 62 if necessary, we have the fol- 
lowing. There is an as-solution x_(t), t G (^2-5^0] to some of the systems (a, ±/) in M 
such that X- (to) = x+(to) = z{so), and X-\(t2-,to] and z\{so—Si, Sq] are homeomorphisms 
onto z{{so — 61, So]) = x_((t2-,^o])) a^-|('^2-,^o) and z\{so — Si, sq) are analytic diffeomor- 
phisms onto z{{so — 5i,So)) = a;_((t2-, ^0))) which is an analytic submanifold which is 
also a semianalytic subset. Moreover, (x„|(t2-, ^o))""'^ ° (-^Ksq — 5i, Sq)) : (sq — 5i,So) 
is an as-curve in M which is an as- diffeomorfism onto its image {t2-,to). We also have 
thatt2- can be chosen such that for eacht^- G {t2-,to], x_ ((t3_,to)) is a neighborhood of 
X-{to) inx^ ((^2-, ^o]) — {^^-(^o)}- The as-curve z{s), s G [so,so+52) satisfies the equation 
a {z{s)) z{s) = i{s)f (z(s)) , s G (sq, So+62) where i(s) > 0, s G (sq, So + 52)- The as-curve 
z{s), s G (sq — 61, So] satisfies the equation a {z{s)) z{s) = i{s)f {z{s)) , s G (so — 61, Sq) 
where i{s) > if x_{t) satisfies {a, f) and t(s) < if x_(t) satisfies {a,—f). Similar 
results hold for as-solutions X-{t), t G (to,i^i]- 

Proof. Part (a) and part (6) are a direct consequence of lemma 14.91 To prove part 
(c) consider the homogeneous IDE {a, 0) with domain R x M and range F where 
a(t,x)(t,x) = a{x)x — tf{x). Consider the solution (t{s), z{s)), s G (so, sq + ^2) to the 
homogeneous system (a,0) where z{s), s G (sq — Si, sq + ^2) is the as-curve considered 
in (6). Since we know from (6) that a {z{s)) i(s) and / {z{s))) are linearly dependent for 
s G (so, So + ^2) we can conclude, using the analyticity of z{s), s G (sq — 5i,So + 62), 
that they must also be linearly dependent for s G (so — Si,so + 62). We can also 
show, by using lemma lOl that there exist 61 such that z {{sq — 61, sq]) — {z{so)} is 
nonempty and locally connected at z{so), more precisely, 2; ((so — ^a, sq)) is a neigh- 
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borhood of z{so) in z{{sq — Si,So\) — {z{so)} for all Sq — S3 G (sq — 5i,So). Moreover, 
Si can be chosen such that z\[so — Si, Sq] '■ [sq — Si, Sq] z ([sq — ^i, Sq]) is an home- 
omorphism, z ((sq — 5i, So)) is an analytic submanifold which is a semianalytic subset 
of M and z : (so — ^i^-So) 2; ((sq — 5i, sq)) is an analytic diffeomorphism. Since 
we have a linear dependence between a and / (^(s)) , s G (sq — 5i, Sq + (^2) 

we have several cases. Assume first that a(2;(s))i(s) = 0, s G (so,So + 52)- Then 
since t(s) 7^ 0, s G (soi^o + ^2) we can conclude that f {z{s)) = 0, s G (sq, Sq + S2)- 
By analyticity of z{s), s G (sq — Si,sq + ^2) we obtain that a {z{s)) z{s) = and 
/ {z{s)) = 0, s G (so — 5i,so + S2). Then to prove (c) in this case we can simply take 
t — s — So + to, t2- — —Si + to and X2-{t) — z{t + so — to)- We can proceed in a simi- 
lar way if we assume that f {z{s)) — 0, s & {so, so + S2)- Let us consider now the case 
where a {z{s)) z{s), s G (sq, Sq + S2) is not identically 0. We are going to show that, after 
conveniently diminishing the value of Si if necessary, there exists a unique as-curve in 
M, A : (so — Si, So) — >■ M such that a {z{s)) z{s) = \{s)f {z{s)) , s E {so — Si,sq). First 
one should take into account that z{s), s G (sq — 5i,so + S2) is analytic then so are 
/ {z{s)) and a {z{s)) z{s), s E {so — Si, So + S2) and therefore they have at most a finite 
number of isolated zeros in a neighborhood of Sq. It is easy to see from the equation 
a {z{s)) z{s) = A(s)/ {z{s)) , s E {sq — Si, Sq + S2) that A(s), or rather its extension for 
complex s, is a meromorphic function in a neighborhood of sq- On the other hand, since 
\{s) = i{s), s G {so,so + ^i) where t{s) is bounded, we have that A(s) cannot have a 
pole at So therefore it must be analytic in a neighborhood of sq- This implies that by 
conveniently diminishing the value of Si and the value of S2 if necessary, we have that 
A(s), s G (sq — 5i, So + '^2) is real analytic. By conveniently diminishing the value of Si if 
necessary we can assume without loss of generality that A(s) 7^ 0, s G (so — 5i, sq). Let us 
assume first that A > 0. Then the result follows by taking t — t{s), s G (so — 5i, sq + S2) 
such that dt/ds — A(s), t{so) — to and t{so — Si) — t2-, which defines t{s), s G (sq — so] 
as an as-curvc in M and also t2-- Since A(s) > 0, s G (sq — 5i,So) we have that t{s) is 
an analytic diffeomorphism from (sq — 5i,So) onto (^2-1^0) which is an as-curve, while 
t : [sq — Si, So] [t2-,^o] is an homeomorphism. Then we can define X2-{t), t G (t2-,^o] 
by X2-{t) = z (s(t)) where s{t) is the inverse of t{s). The rest of the proof can be per- 
formed in a similar way. ■ 

The previous theorem says, in particular, that an as-solution x^(t), t G [to,ti) to a 
given IDE (a,/) gives rise to an as-solution X-{t), t G (t2-,^o] of (a,/) or (a,—/) such 
that x+(to) = X-{to) and graph a;_ U graph a;+ = graphs for some as-curve z{s), s G 
(sq — Si, So + S2), which is a solution to the homogeneous system a{z)z — if{z), s G 
(sq— 5i, So)U(so, S0+S2), and x_, x+ are reparametrizations of z\{so—Si, Sq], 2;|[so, So+52)- 

For a given IDE (a,/) with domain M and range F, an as-solution x{t), t G [to,ti) 
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(t G (to)^i]) t G (to,^i), t £ [^0)^i]); to (a, ±/), in M is, by definition, an as-solution 
to some of tlie systems (a, /) or (a, — /), in M. Tlie notion of a pas-solution or an Ics- 
solution to (a, ±/), in M is defined by tlie condition that each as-piece is an as-solution 
to (a, ±/) in M. 

We have the following extension theorem 

Theorem 4.17 Let (a, /) be an IDE with domain M and range F. Let x{t), t G [toj^i] 
he a continuous pas-solution to (a, /), so its graph is compact. Then x{t) can be extended 
to a continuous pas-solution x{t), t G [to)^i]) (a, ±/) for some < t^ < ti < ti. 

Proof. If X is a constant the proof is immediate. Let us assume that x is not a 
constant. We have an as-piece of x of the type x|[to,to + S2). As a consequence of 
theorem I4.1fjl we have that x|[to,to + ^2) can be extended to a continuous pas-solution, 
say u{t), t G [t2-,to -\- S2), to (a, ±/), where u+ = x| [to, to -\- 62) and U- : {[t2-,to] M 
are the as-pieces. We can glue the as-piece to x and take to = ^2- to obtain a 
continuous pas-solution to (a, ±/) defined in [io,ti) which extends x. We can proceed 
now in a similar way to obtain, in turn, an extension of this solution to a pas-continuous 
solution to (a, ±/) defined in some interval [to,^i]- ^ 

Corollary 4.18 Let (a, 0) be an homogeneous IDE with domain M and range F. Let 
x{t), t G [to,^i] be a continuous pas-solution to {a,f), so its graph is compact. Then 
x{t) can be extended to a continuous pas-solution x{t), t G [^05^1], to (a, 0) for some 

to<to<ti < ti. 

5 The example of the symmetric elastic sphere 

As we have said before the main purpose of this paper is not to pursue the investigation 
of how to obtain systematically desingularizations of given IDE by using systematic 
procedures to desingularize given closed analytic sets and then apply all this to solve a 
given IDE by a systematic procedure, but rather to show how having a desingularization 
of a given IDE may help to understand its solutions. In the next section we study 
an example from nonholonomic mechanics, namely the symmetric elastic sphere, and 
show using a desingularization of the equations of motion reduced by the symmetry, 
which constitute an IDE, how to study the dynamics and solve completely the system 
by quadratures. 
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Preliminaries. Continuing with tlie description of the symmetric elastic sphere given 
in section we shall now give a more precise description of this system as a nonholo- 
nomic system. For references on nonholonomic mechanics related to rolling spheres see 
[3 El 113 EH EH] and references therein. 

A rigid sphere rolling on a plane can be modelled as a nonholonomic system on the 
group 5*0(3) X where, for a given element {A, x) G 5*0(3) x R^, A represents a rigid 
rotation and x the position of the point of contact of the sphere with the plane. Now 
we will assume that 5*0(3) x is also a good configuration space for an elastic sphere 
where deformations are small and concentrated only near the area of contact which is a 
small circle whose center has a position given by x. The kinematics of this system can be 
conveniently described as follows. We assume that there is an ortonormal system fixed 
in the space, say (61,62,63), 61 = (1,0,0), 62 = (0,1,0), 63 = (0,0,1), then we have a 
basis moving with the body, {Aei, Ae2, Ae^), where A = A{t). We introduce the variable 
z E S^, given hy z = Ae^. The spatial angular velocity u can be written u = vqz + z x z, 
so Vq = {uj, z) is the component of 00 along z. The nonholonomic constraint is given 
by the nonsliding condition. Since the area of contact of the sphere with the plane 
is a circle we must add to the usual nonsliding condition u x r63 = x for the rigid 
rolling sphere (see, for instance, where Lagrange-D'Alembert-Poincare equations 

for the rigid sphere have been written) where r is the radius of the sphere, the extra 
condition that the vertical component of the spatial angular velocity is 0, that is, = 0. 

We are going to assume that the center of mass coincides with the center of the sphere 
and that the principal axis of inertia are [Aei, Ae2, Ae^). The three principal moments 
of inertia of the sphere are Ji, I2, I3, and we are going to assume that Ji = I2. We 
introduce the adimensional quantities a = /3//1 and (3 = Mr^//i, where M is the mass 
of the sphere. The Lagrangian of the system is given by the kinetic energy. 



where x is the velocity of the center of the sphere. The nonholonomic constraint is given 
by i = ti; X r63 and ^3 = 0, and using this we can conclude that the kinetic energy of 
the actual motion of the symmetric sphere is given by 



As we have said in the Introduction the addition of the extra condition = introduces 
an extra singularity in the reduced system, which is an IDE, and we need to apply our 



hz' + Ihv', + Imx' 
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desingularization procedure to obtain a single equivalent differential equation describing 
the system, in the sense of theorems 14 . 1 41 and 14 . 1 51 We will show that the desingularizing 
manifold containing the essential dynamics, in which this differential equation is defined, 
is diffeomorphic to S"^ x S^. Another interesting feature of the desingularization method 
for the example of the symmetric elastic sphere is that integrability by quadratures ap- 
pears in a natural way, in terms of angular coordinates of S"^ x S^. 



The IDE for the symmetric elastic sphere. As a result of reduction by the sym- 
metry techniques, in this case reduction by the subgroup 5*0(2) x we obtain the 
following system of Lagrange-D'Alembert-Poincare equations, which is an IDE, 

{a + (3){zx e^)vQ + [1+ < z,e^> Vii- 

{a + [3)vq < z,e^> {z X z) = (5.1) 
Vq < z,e3 > + < z X z,e3 > = 0. (5.2) 

Here V represents the Levi-Civita connection on S*^ with respect to the standard met- 
ric. This is a consequence of the methods developed in after some more or less 
straightforward calculations which we will not explain here. The previous Lagrange- 
D'Alembert-Poincare equations are derived under the assumption ^3 7^ because the so 
called dimension assumption adopted in [S] is not satisfied for the whole manifold S"^. 
Nevertheless, by continuity, the equations ()5.H) . ()5.2|) are also satisfied by the motion of 
the rolling ball for z^ = 0. Without using the derivation of the Lagrange-D'Alembert- 
Poincare equations, the careful reader may want to check directly that the previous 
system of equations, or equivalently, the system of equations f|5.8|) - ()5.15p . is equivalent 
to balance of momentum plus the condition uj^ = 0. 

Since Vii = z x [z x z) hj taking the inner product of (j5.H) with z and using (j5.2j) we 
get, at least for < 2;, 63 > 7^ 0, 

= ±{{l + P)z' + {a + P)vl), (5.3) 

from which one deduces 

2e = (1 + I3)z^ + (a + (3)vl (5.4) 

where e represents the normalized energy. This equation represents conservation of en- 
ergy, as one can check more directly by looking at the expression of the kinetic energy 
E given at the beginning of this section. We shall assume from now on that e > 0, 
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otherwise the motion is triviaL 



We have the following equations to be satisfied for the symmetric elastic sphere in vari- 
ables {z, u) where z = v and v x z = u, so the variable vq does not appears, 

(1 + /5) < 2:, 63 >< -it, 63 X z > +(a + /5) < M, 63 >^ = (5.5) 
(1 + ^) < z,e3 >^ + (a + < u,e3 >^ -2e < 2,63 >^ = 0, (5.6) 

Equation (j5.5|) is obtained by taking the inner product of (j5.ip with 63 and using (j5.2|) 
while equation ()5.6p is obtained from equation ()5.4|) and equation ()5.2p . 

Equations ()5.5|) and ()5.6|) involve the variables (2;, u) G TS"^ and we have a natural 
inclusion TS"^ C S*^ x R^, where 

TS^ = {(z, m) G ^2 X .^z^u >= 0}. 

Equations (j5.5p . (|5.6|) form an IDE in the manifold TS"^. By adding the equation z"^ = 1, 
we obtain an equivalent IDE in the variables {z, m) G x M^. Let us include, in addition, 
the equation of conservation of energy ()5.4|) . written in terms of the variables {z,u,vo), 
as follows 

2e = (1 + + (a + I3)vl (5.7) 

In other words, we are going to study the system of equations given by (|5.5|) . (j5.6|) . 
(j5.7|) for a fixed e > 0. Of course equations (j5.6p and (|5.7p taking into account (j5.2p 
are redundant for Z3 7^ 0, but for z^ = the system given by the equations ()5.5|) . ()5.6|) 
includes solutions of the type z(t) = const, with z^ = 0, where fo takes any given value 
and, since u = for this kind of motion, the energy is given by 2e = (a + /5)(fo)^ and 
therefore the condition that the energy must have a fixed value will not be satisfied. Of 
course we can study with our methods both, the system given by (j5.5|) . (j5.6j) . (j5.7p and 
also the system given by (j5.5p . (j5.6|) . but we will chose to study just the first of them, 
for simplicity. 

The IDE for the symmetric elastic sphere in the standard form. Considering 
that z^ = 1 and z = z x u, we must have < z,u >= 0. Using what was said in the 
previous paragraph, and according to theorem 12.11 we can write the system of equations 
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for the symmetric elastic sphere in variables {z, u, Vq) G 



X 



as follows, 





= Z2U3 - Z3U2 




Z2 


= Z3UI - ZIU3 




Z3 


= ZIU2 - Z2UI 







= (1+ P)Z3{-Z2U1 + Z1U2) + 


(a + P)ul 





= il+(3)zl{ul + ul + ul) + { 


a + I3)ul- 2ezl 





= zf + zl + zl-l 







= ZiUi + Z2U2 + -23^3 







= 2e- (a + /3)wo' 





The system ()5.8|) - ()5.15|) can be written in the form 

a{X)X = f{X), 

with X = {z,u,vo), where, 



a{z,u,Vo) 



1 





1 



1 




























f{z,u,vo) 











Z2U3 - Z3U2 
Z3UI - ZIU3 
ZIU2 - Z2UI 
~{a + l3)ul 
(1 + P)zl{u\ + ul + ul) + (a + P)ul 
zl + zl + zl-l 

ZiUi + Z2U2 + -23^3 









il+P)z2Z3 {l+^)ziZ3 













2ezl 



2e 



(5.8) 
(5.9) 
(5.10) 
(5.11) 
(5.12) 
(5.13) 
(5.14) 
(5.15) 



Application of the algorithm. We will work on the manifold M = M^, where 
{zi, Z2, Z3,ui,U2,U3,vo) G M7 are independent variables. Then our IDE is given by 
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equations We can easily see that kr = 4, 5*4 (M) = M, L/^{M) = Mq, 

Ml = M - L4(M), M2 = 0. Now we shall describe Mq by equations. Let 



V^l = -{1+ P)Z2Z3 
V92 = (1 + P)ZiZ3 

i/i = (1 + P)zl{ul + ul + ul) + {a + P)ul - 2ezl 

1^2 = zl + zj + zl - 1 

U3 = ziui + Z2U2 + 2:3^3 



(5.16) 
(5.17) 
(5.18) 
(5.19) 
(5.20) 
(5.21) 



z/4 = 2e - (1 + P)u^ - (a + (3)vl. 



As we know Mq = Li{M) is given by the condition that rank[a, /] < 4. Let 



Moa = {y^l = 0,^2= 0} 

= {Z3 = 0}U{ZI = 0,Z2=0} 

Mob = {ui = 0,1^2 = 0, z/3 = 0, z/4 = 0}, . 



(5.22) 
(5.23) 
(5.24) 



Then we can easily see that Mq = Moa U Mof,. The desingularization M-^ of Mq will be 
the disjoint union of the desingularizations of Mqo and Mq^. 

The desingularization of Moa can be described by = {z^ = 0} \J{zi = 0,Z2 = 0}, 
where \_\ means disjoint union and the projection ttq is the identity on each disjoint piece 
of M^. One can see using ()5.8p - ()5.15|) that the lifted system (ai, fi)\{z3 = 0} satisfies 
Z3 = 0, ^3 = 0, zf + Z2 = 1, which implies z = 0, and also, since u = z x z, that 
u = 0. This describes the motion completely. It consists of the rolling of the sphere 
with z{t) = {ziQ, Z2Q, 0) fixed and the z component of the angular velocity vq satisfies 
2e = (a + P){vqY. The lifted system (ai,/i)|{zi = 0,Z2 = 0} satisfies zi = 0, Z2 = 0, 
Z3 = ±1, therefore z = 0, and then u = 0, which contradicts equation z/i = 0, because 
we have assumed e > 0. So there is no motion, that is, no solution, for the system 



Now we will desingularize Mob- We are going to see that Mob is in fact a nonsingular 
manifold. More precisely, we will define the desingularizing manifold M^ by equations in 
the variables {z, u, Vo), with voz^ = U3, from ()5.2p . For simplicity, we call fj, = 2e/ (1+/3) > 
and A = {a + (3) / {1 + (3) > 0, from now on. Then we have the following equations 
defining the nonsingular manifold M^, 



(ai,/i)|{^i = 0,Z2 = 0}. 
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= - vqZs (5.25) 

= ul + ul + ul + Xvl - (5.26) 

= zl + z^ + zl-1 (5.27) 

= ziUi + Z2U2 + Zsu-i- (5.28) 

See Appendix A for a proof that the system above defines a nonsingular manifold. The 

map TTo : — * M is then given by the restriction of the identity {z, u, vq) —>■ {z, u, vq) 
to and one can check that the image of ttq is precisely M^h- 

According to theorem 12.11 the system lifted to has the same solutions as the system 
given by the equations 

zi = Z2U3 - Z3U2 (5.29) 

Z2 = Z3U1 - zius (5.30) 

Z3 = Z1U2 - Z2U1 (5.31) 

2:2^1 — Z1U2 = XvqUs (5.32) 

= M3 — V0Z3 (5.33) 

= ul + ul + uj + Xvl- jj (5.34) 

= zl + zl + zl-1 (5.35) 

= ziui + Z2U2 + z-iU-i. (5.36) 



More precisely, if we define 



a{z,u,VQ) 



10 
10 



1 















Z2 

























-zi 













we see that the system ()5.29|) - ()5.36|) is in the form a{y)y 
it is an IDE in standard form with domain M7 and range H 



Z2U3 - Z3U2 
Z3U1 - Z1U3 
Z1U2 - Z2U1 

\v0U3 
U3 - V0Z3 
u'l + ul + ul + Xvq — fi 

zf + Z2+zl — l 
ZiUi + Z2U2 + Z3U3 



-- f{y), with y = {z, u, Vq), so 
\ and our IDE in is given 
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by the restriction (ai, /i) = (a, /)|M^^. In order to continue with the algorithm, we shall 
find explicitly the lifted system (ai,/i). By diff'erentiating the equations (|5.25|) - (j5.28|) . 
eliminating the redundant equation ziZi + Z2Z2 + z^z^ = 0, also realizing appropriate 
hnear operations in the range space, that is, using projections as explained in section |2l 
and also according to theorem 12.11 we have the following system with domain MJ and 
range M}^, which is also equivalent to our system (ai, /i). 





= Z2U3 - Z3U2 


(5.37) 




= Z3U1 - Z1U3 


(5.38) 




= Z1U2 - Z2U1 


(5.39) 


Z2U1 - Z1U2 




(5.40) 


Ziili + Z2U2 + -23^3 


= 


(5.41) 


uiiii + U2U2 + + XvqVq 


= 


(5.42) 


U3 - Z3V0 


= V0ZIU2 - V0Z2UI 


(5.43) 





= U3- V0Z3 


(5.44) 





= ul + ul + ul + Xvl — fj, 


(5.45) 





= zl + zl + zl-1 


(5.46) 





= ZiUi + Z2U2 + Z3U3. 


(5.47) 



This system is still not completely desingularized. One can check by direct calculation 
that it can be desingularized in two more iterations of the algorithm. However, in this 
example there is an interesting alternative to find the solutions, which starts with a 
precise description of the manifold Ml. We prefer this alternative because having an 
identification of Ml also helps to understand the dynamics in a direct way, as we will 
see soon. This is precisely one of the points of the present work, namely, to show how 
in examples the desingularization method may helps to understand the dynamics of a 
given system. 

Identification of Ml. The manifold Ml is given by the equations ()5.25|) - ()5.28|) in 

the space of the variables {zi, Z2, ^3, Mi, U2, M3, t'o), as we have seen before. The equation 
(j5.28j) tells us that u is a vector tangent to the 2-sphere S'^, given by 2;^ — 1 = 0. 
Heuristically, for each 2; e we consider the 3 dimensional space T^S"^ x where 
represents a line normal to the sphere at z G 5*2, so the collection of all Rz is a trivial real 
line vector bundle with base S'^. Equation ()5.25|) is a plane containing the origin = 0^ 
since Z3 is fixed once z is fixed. Equation ()5.26|) gives an ellipsoid. The intersection 
of the plane with the ellipsoid is an ellipse. Therefore Ml must be some fiber bundle 
with fiber S"^ and base S'^. Using all this and some imagination we can see that it is. 
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in fact, the trivial bundle S"^ x S^, moreover, we have the following parametrization of 
in variables {6, ip, ip). In any case, this assertion can be easily checked after some 



straightforward calculations. 

Zl = sin6'cosv9 (5.48) 

Z2 = sin 6^ sin (5.49) 

Z3 = cose (5.50) 

Ml = —acos{(f — ip)cos^9cos(f — bsm(Lp — -ip)smLp (5.51) 

M2 = —a cos{ip — ip) cos^ 9 sin {p + b sm{{p — ip) cos ip (5.52) 

Ms = acos((y9 — ■?/;) cos6'sin6' (5.53) 

Mo = a cos(v9 — t/;) sin 6^, (5.54) 



where 

Q = \/ , ■ 2/1 T^' b=^'. 

V A sm 6 + cos^ 9 

In other words, by some straightforward calculations we can check that {zi, Z2, Z3, Mi, M2, M3, Vq) 
in coordinates {9, (p, ip) satisfies ()5.25|) - ()5.28|) . 

We can see that equations ()5.48|1 - ()5.54|1 define a diffeomorphism f : S"^ x 
f {z , {cos ip , sin ip)) = {z,u,vo), which gives the desired identification of M^. See Ap- 
pendix B for a proof that / is a diffeomorphism. 



The differential equation for the symmetric elastic sphere in in variables 

{9,(p,ip). Considering the parametrization for given by ()5.48|) - ()5.54j) we get the 
equations for ()5.29p - ()5.36|) in coordinates {9,ip,ip), 

cos 9 cos ip9 — sin 9 sin (fip = a cos 9 sin ip cos{ip — ip) — (5.55) 

—b cos 9 cos ip sin{{p — ip) 

cos 9 sin ip9 + sin 9 cos (pip = —a cos 9 cos ip cos{(p — ip) — (5.56) 

—6 cos 9 sin ip sin{(p — ip) 

— sin99 = b sin 9 sin{ip — tp) (5.57) 

a sin cos^ 6' cos(v9 — ?/')0 — (5.58) 
—b sin 9 cos{lp — %Ij){lp — Ip) = Xa^ cos^ (ip — ip) sin^ 9 cos 9 
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If sin 6^ 7^ the system (|5.55|) - (j5.58|) becomes 



e = -bsm{ip-ip) (5.59) 
cos 6* 

ip = —a -cos{ip — yj) (5.60) 

sin 9 

ip = a cos{ip - ip)^^?^ ( - - l) , (5.61) 
sm V a / 



or equivalently, 



9 = -bsm{^-ij) (5.62) 

(f = —a cos(v9 — ip) (5.63) 

sin 9 

cos 9 

Ip = {b — a) -cos{(p — ip) (5.64) 

sine* 

It can be easily seen that this system can be integrated by quadratures. For instance, if 
we call w = ip — Ip, we can write (|5.62|) - (j5.64j) as a planar system in coordinates {9, w), 

9 = —bsinw (5.65) 
cos 9 

ill = —b -cosw, (5.66) 

sm 9 

which in turn leads to the separable equation 

d9 ^ .X 

-— = tant'tanw. (5.67) 

aw 

Of course the system ()5.55j) - ()5.58|1 is still an analytic IDE on the analytic manifold 
5^ X which is of constant rank for sin 6' 7^ 0, and the rank changes for sin 6' = 0. 
So we should continue the desingularization process, which is not difficult but we are 
not going to explain the details. Instead, we simply observe that it is not difficult to 
see, alternatively, that the only solution with some initial condition compatible with 
the system and involving the condition sin 6' = 0, that is, initial condition of the type 
(^10) Z30, Mio, U20, M30) ^00) = (0, 0, ±1, Mio, U20, 0, 0), consists of a uniform circular mo- 
tion of 2 on a vertical plane perpendicular to the constant vector {ui(t) , U2{t) , u^it)) = 
(mio, M20, 0), while fo(t) = 0. This is also consistent with physical reasoning. 



We are not going to perform a detailed study of the dynamics of the symmetric elastic 
sphere in this paper, although it is clear that having the kind of explicit solutions that 
we have found clearly helps to study the standard dynamical systems questions about 
equilibria, linear and nonlinear stability, bifurcations, and, also control of the system, 
with potential applications to robotics. A future work in this direction is being planed. 
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Appendix A. We will prove that is a nonsingular 3-dimensional manifold. We 
must prove that the matrix 









-Vo 








1 


-Z3 











Ui 








Zl 


Z2 
















Ui 


U2 


U3 


Zl 


Z2 


Z3 






(5.68) 



has rank 4, for all {z, u, Vq) G Ml. 



1. If zz = ±1, zi^ Z'2.^U3^vo^ 0, (iti, U2) 7^ 0, we get 

■ 1 ±1 

W2 

^1 

ui U2 ^ ^1 

which has rank 4. 



2. If ^3 = 0, M3 = 0, (2:1, Z2) 7^ 0, + Z2U2 = 0, we have 



-?7o 1 

Ui U2 Xvq 

Zl Z2 

Ui U2 ^1 ^2 



(a) If = 0, then U2) 7^ 0. For instance, if 7^ we have the subdeterminant 

-ui S S ^0. 



In fact, if 



then 



1 

Ml 

^2 

ui U2 Zl 



Zl Z2 
Ui U2 



Zl Z2 
Ui U2 



0, 



-Z2U1 
ZlUi 



ZIU2 
Z2U2 
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and since zf + Z2 — 1, we have {ui, U2) = 0, a contradiction. A similar result 
holds if 7^ 0. 

(b) If vq ^ 0, and, say z\ ^ 0, we have 



-vo 


zi 
ui 





Ml 














Zl 
Ui 



-Vo 

Zl 



-Xvlzf + 0. 



A similar result holds if 2^2 7^ 0. 



3. Assume now that < \z-i\ < 1; (-21,-22) 7^ 0- 



(a) If Vq = 0; 1*3 = 0, {ui, U2) 7^ 0, ZiUi + -^21*2 = 0, we 



have 


















1 


-zs 











Ui 










Zl 


Z2 


Z3 














Ui 


U2 





^1 


Z2 


Z3 






If, for instance, Ui ^ 0, we have 












-^3 








Ui 





Zl 


Z2 








Ui 


U2 


^1 






Zl Z2 
Ui U2 



In fact, if 



Zl Z2 
Ui U2 







then — -2;2tii+-2iii2 = and this together with Z1U1+Z2U2 
or (-^1,-22) = 0, a contradiction. Similarly if ^2 7^ 0. 



gives {ui, U2) — 
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(b) If Vq 7^ 0; M3 7^ 0, ZiUi + Z2U2 7^ 0, {ui,U2) 7^ 0; in this case the matrix 
(|5.68|) also has rank 4. In fact, (2:1,^25-23) and {ui^U2,u^) must be linearly- 
independent, otherwise the condition ziUi + Z2U2 + ^^s^s = would imply 
M = 0. Then the last two rows of ()5.68|) are linearly independent. Using 
this we can see easily that the second row cannot be a linear combination of 
the last two rows, so the last three rows are linearly independent. Let 
i = 1, 2, 3, 4 represent the four rows. If the rank of the matrix is not 4 then 
we must have -Ri = X-R4 + i'Rs + V^-R2- This leads to the systems of equations 

XUi + i>zi = 
XU2 + ipZ2 = 

XU3 + i)Z-i = ~Vo 

and 

XZi + (fui = 

XZ2 + ^U2 = 
XZ3 + ifU^ = 1. 

Using < z,u >= we can deduce from the first system that = —VqU^ and 
since ^3 = VqZs we obtain xu^ = —^qZ^. On the other hand, using < z,u >= 
we can deduce from the second system that x = ^3. We can conclude that 

= —Vq, which implies u = vq = 0, a. contradiction. 
Therefore, also in this case the rank is 4. 

We can conclude that is a nonsingular manifold of dimension 3, defined regularly in 
by the equations (jE2il)-(EI2Hl)- 

Appendix B. The proof that / is a diffeomorphism goes as follows. 
1. Injectivity 
(a) ^3 7^ ±1. 

For given {z, {cosip, sinip)) G S*^ x S^, we have {z, u, Vq) = f{z, {cosip, sin ■?/')). 
We have that (6*, y?) G (0, vr) x [0, 2it) is uniquely determined, moreover, < 
sin 6' and — 1 < cos^ < 1 are determined. Therefore from equation (|5.54|) 
cos{ip — ilj) is determined. Using ()5.51|) and ()5.52|) we obtain sm{ip — ip). 
Since z determines (cosy?, sin yj), we get ip G [0,27r) uniquely determined. 
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(b) = ±1. 

In this case sin^ 
obtain 



0, cos6l = ±1, a = 6 = y/I. Using (jK3T]) and we 



Ml 
U2 



—y/Jl cos ip 



which determines ip G [0, 2tt). 



2. We must prove now that the rank of the tangent map T(^z, (cos i(>,sim(>))f is 3 for 
all {z, {cosipjSinip)) G S*^ x S^. We have that the Jacobian matrix calculated at 
(6*, (f, tp) of the map given by ()5.48|l - ()5.54j) is 



a' cos(v? 



where we have written 



cos cos if 

cos 6 sin if 
— sin 9 
J41 

ip) sin ^ + a cos(v9 — ip) cos 6' 



— sin 6 sin 
sin 9 cos 


Jg2 









-a sin(<^ — '?/') sin ^ a sin (y) — ■?/') sin ^ 



J41 

<^42 

Jei 



—a' cos((y9 — ip) cos^ COS + 2a cos((y9 — ip) cos 6* sin 6* cos ip 

a sin((y5 — ip) cos^ ^ cos + a cos(y9 — ip) cos^ ^ sin Lp— 

—b cos{(f — ip) sin Lp — h si\i{(p — ip) cos ip 

—a sin(y9 — lp) cos^ ^ cos ip + h cos{ip — ^Z") sin ip 

—a' cos(v9 — '?/') cos^ Osuiip + 2a cos{ip — ip) cos 9 sin 6^ sin ip 

a sin((y9 — lp) cos^ 6' sin — a cos(v9 — ip) cos^ 6' cos 93+ 

+b cos(v9 — ^Z)) cos ip — b sm{ip — ip) sin ip 

—a sm{ip — lp) cos^ 6 sin yj— 

—6 cos((y9 — '?/') cos V9 

a' cos(v9 — '?/') cos 6* sin 6* + a cos(v9 — ip) [— sin^ 9 + cos^ 6*] 
—a sin(y9 — t/;) cos 9 sin 6' 
a sin(y) — lp) cos sin 9 



with a' = da/d9. 



(a) If 2:3 = ±1 then ^ = or = vr. Let us consider the case 9 = only since 
the case 9 = n can be treated in a similar way. We cannot use the polar co- 
ordinates representations, in fact, we can see that the rank of J{0,ip,ip) and 
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J{tt, ip) is less than 3 for all (</?, ■0). On the other hand, we can easily check 
that for fixed ip, f{9(t), ip, ip) represent the same point in the manifold for 
6(t) = and arbitrary, which is the point {0,0,1, —bcosip, —bsinip, 0,0). 
Let, for instance, 9{t) = t, then the tangent vector to the curve f{9{t),ip,ip) 
att — O, is (cos 99, sinip, 0, 0, 0, bcos{(f — bcos{ip — ip)) ■ For ip — and (f — 
7r/2 we obtain the hnearly independent vectors. A = (1, 0, 0, 0, 0, b cos -0, b cos 0) 
and B = (0, 1, 0, 0, 0, 6sin-0, bsini/j) , respectively. On the other hand, the tan- 
gent vector to the curve f{0,0,ilj + t) at t = 0, which is the last column of 
J(0, 0, ip), is C = (0, 0, 0, bsinip, —bcosip, 0, 0). The three vectors A, B, C are 
linearly independent for all -0, which shows that the rank of T(^e3,{cosi),sinip))f 
is 3 for all ^0. 

(b) If zz — 0, cos^ = 0, then 9 — 7r/2 and a — -s/n/X- We get in this case 



J(7r/2,(^,V') = 






— sin (/? 










COS(^ 







-1 













— 6sin(2(y9 - 


-0) 


bcos{(p> — ip) sin 93 





bcos{2(p — 


0) 


—b cos((^ — V^) cos ip 


a cos(i^ — 0) 













— asin((/7 — 


0) 


asin((/7 — 0) 



i. If sin 9? = 0, then cos 9? = ±1. 

Case I: sin((/7 — ■0) ^ 0, then the rank of J is 3. 

Case II: sin(99 — 0) = 0, then cos((/7— ■0) = ±1 and — 6cos(</7— ■0) cost/? ^ 0; 
then the rank of J is 3. 

ii. If sin(/9 = ±1, then cos(/9 = 0. 

Case I: sin(9? — 0) 7^ 0, then the rank of J is 3. 

Case II: sin((^ — 0^) = 0, then cos(9? — 0) — ±1 and — 6cos((^— 0^) sin 99 7^ 0; 
then the rank of J is 3. 

iii. If < I siny?| < 1, 

Case I: sin((y9 — -0) 7^ 0, then the rank of J is 3. 

Case II: sin(v9 — -0) =0, then cos{ip — ■0) = ±1 and b cos((^ — ip) sin 7^ 0; 
then the rank of J is 3. 



(c) If < |2;3| < 1, then we have 

20S 9 cos 09 — sin ^ sin 09 /,./,/ ^ 

„ . . „ =cos^sm^7tO; 

cos 9 sm (y9 sm 9 cos 

then the rank of J will be 3 if Jjs 7^ 0, for at least one i, 4 < i < 7. 



■hi 
J21 



J12 
J 22 
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i. If J73 7^ 0, then the rank of J is 3. 



ii. If J73 = 0, then sm{ip — ip) = and cos((/9 — ip) — ±1. Moreover, J43 = 
±6 sin and J53 = ±a cos^ 9 sin ip ±b cos 
Case I: If sin 7^ 0, then the rank of J is 3. 

Case II: If sin (p — 0, cos cp — ±1 and J53 = ±6 7^ 0, then the rank of J is 
3. 
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